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APPENDIX B 

CLUSTER analysis OF SKYLAB RADSCAT DATA 
by A, Sobti, Satnut Siriburi, and R. K, Moore 

ABSTRACT 

The Skylab S-IPS radiomei:er/scaf•^erometer produced ferrain responses with 

2 

various polarizations and observation angles for cells of 100 to 400 km area. 
Classification of the observations into natural categories (for the sensor) was achieved 
by K-means and spatial clustering algorithms. The principal land use in each 
observation cell, or in representative- sample cells, was determined by use of maps 
and Skylab S-190 photographyo Little correlation was found between the clusters 
natural to the data and the land-use categories determined in this way, although water 
and land surfaces were usually separable and some gross trends could be inferred. 
Apparently natural clusters in microwave data are sufficiently affected by other factors 
to render such gross- resolution data ineffective for land use modelling. ‘ 
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INTRODUCTION 


The Skylab radiometer/scatterometer operating at 13.9 GHz had a resolution 

cell ranging from approximately 100 sq, kms, at nadir incidence to approximately 400 

sq, km, at 50° incidence. To determine if with such gross resolution the microwave 

response could be partitioned to identify or discriminate between major terrain features 

as discerned by maps and imagery, some clustering exercises were conducted. The 

Skylab radiometer/scatterometer provided backscattering coefficient and radiometric 

temperature at various angles of incidence and for various polarizations. The in~track 

modes provided multiple”angle observations of a- target and the cross~track contiguous 

modes provided either a pair of radiometer or scatterometer measurements or one of 

each. The clustering exercises reported here are from data obtained with the in~track- 

contiguous and cross~track- contiguous pitch-offset 29° modes. The number of data 

points included in each clustering exercise varied and in one case exceeded 7000 points, 

1 

The principal procedure used was K-means clustering . The computations 

were performed on a Honeywell 635 computer. Another clustering procedure, called 

2 

spatial clustering was attempted for one set of data. This latter procedure has been 
used quite effectively in the clustering of images from ERTS, The spatial clustering was 
performed on a PDF 15/20 computer. There was no significant difference in the ability 
to partition the data space into clusters corresponding to major terrain types. 
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MacQueen, J., "Some Methods for Classification and Analysis of Multivariate 
Observations," 5th Berkeley Symposium on Probability and Statistics, 1972. 

Haralick, R. M. and I, Dinstein, "A Spatial Clustering Procedure for Multi-Image 
Data," IEEE Trans. , v. CAS-22, May, 1975. 
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The Skylab Radiometer/Scatferometer Experimeni- 


Skylab had on board a radiometer and scatterometer operating at 13,9 GHz 
which shared the same antenna and much of the rf hardware. The two sensors could 
be operated in near simultaneity; the sequence of observations and the various sensor 
configurations were predesigned into four possible modes which were selectable by 
the astronaut. The operation and characteristics of these sensors has already been 
described in numerous NASA publications, journals and technical reports, so only 
the portions germaine to this study will be reported. 

Of the four modes possible, data from only two are reported in this study. 

The non- contiguous modes were used primarily over the ocean and data from these 
were not considering in the clustering exercises. The in-track“contiguous mode pro- 
vided backscattering coefficient data at approximately 46°, 40°, 33°, 17°, and 3° 
off nadir. The targets with observations at all five angles of incidence were separated 
by approximately 25 kilometers and lay in the along-track direction. Radiometer 
measurements were made at angles from 46° to 3° where scatterometer measurements 
were not being made. These measurements were, however, not included' in our data 
base . The cross-track-contiguous mode provided 1 2 measurements between + 1 1 ° of 
a scan initial point. The scan of the antenna was in the across-track direction and the 
scan could be initialized at 0°, 15°, 30*’ *, or 40° in the pitch or roll direction. This 
cross- track- contiguous (CTC) mode could provide either a combination of a radiometer 
and a scatterometer measurement or two of either sensor with orthogonal polarizations. 

The pitch-off-set 30° submode of the CTC provided 12 measurement pairs with a very 
small change in incidence angles across the scan. This submode was the only one selected 
because the number of data points in the particular angle range (30°~33°) was largest 
for conducting the clustering exercises. 

The targets viewed include areas in North America, South America and ocean 
surfaces. The South American targets were predominantly the tropical rain forest and 
savannah regions in Brazil. Maps and spacecraft imagery were used to determine the 
terrain composition of each cluster. Ideally a cluster based upon microwave responses 
should have been from only a singular terrain category for perfect discrimination or 
clustering. 
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Clustering Procedures 


Two clustering procedures were attempted. The K~means clustering algorithm 
was used on ail the data sets reported in this study. Spatial clustering was used as a 
comparison for one long ITC pass from Texas to Maine on day 253. The two methods 
of clustering are thoroughly described by references 1 for K-means and 2 for spatial 
clustering. The philosophies of the two clustering procedures are distinct and con- 
sequently the clusters formed are not the same, 

K“Means Clustering 

This procedure attempts to break up the data space into similarity groups. 

The procedure is not intended to find some unique, definitive grouping, but rather 
simply to aid the investigator in obtaining a qualitative and quantitative under- 
standing of large amounts of N-dimensional data. The data base included for this 
study consisted of (1) all CTC pitch 29° data where the clustering was performed on 
the radiometer (reduced to emissivity) and scatterometer data for vertical polariza- 
tion; (2) ITC data using all five angles for scatterometer measurements (3) ITC data 
using the middle three angles for scatterometer measurements, and (4) CTC pitch 29° 
data (scatterometer only). 

The procedure involves two parameters; C for "coarsening" and R for "refine- 
ment" . The data base is first converted so that all the data lie within a unit n 
dimensional sphere. This is a scaling procedure and the data can then be weighted by 
assigning weights to individual variables. The program starts with a user-specified 
value of K, and takes the first K points in the sample as the initial means. This selection 
can, however, be forced to select any particular points as the initial means (or seeds, 
as-they are sometimes called). The clustering procedure is an iterative one where the 
nearest sample (in measurement space) based upon, a distance criterion (usually square 
of the difference in magnitudes) is compared to the threshold values (C) selected; if 
the proximity is close enough, the sample Is included in the group and a new mean 
computed. If the sample is greater than a distance R from each cluster, it Is allowed to 
be another seed. In doing this, the group means are updated and each point iteratively 
assigned to groups or allowed to form its own cluster which can attract other samples. 

The final classification is performed with the means of the clusters obtained by the above 
iteration process. This classification is again based upon of least-distance criterion. 
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Figure B-1 .The Non~Supervised Spatial Processing Approach 







Non” Supervised SpaHal Clustering 

The spatial clustering algorithm differs from the usual clustering algorithms 
in that contextual information is used in the form of the spatial location of target 
cells. This algorithm was developed for analysis of image-like data where the 
partitioning of data into like regions must take into account the spatial distribution 
of each cell. The algorithm is carried out in two phases. 

In the first phase, the spatial information of the image-like data base is 
used to determine regions of homogeneity. To create an image-like data base from 
data collected by the in-track modes of the S-193 sensors, the target ceil was re- 
peated in one direction (orthogonal to the track) so that there was a grid produced 
which had all columns equal. The clustering then only affected the' rows. The regions 
of homogeneity are found by a boundary enhancement technique. This is shown in 
Figure B-1 as the gradient operator. The gradient image so produced is then processed, 
through a threshold detector and the output image should then contain regions of 
homogeneity separated by boundaries where the gradient exceeded the threshold. 

The next step is labelling each homogeneous region. 

The second phase of the algorithm consists of an iterative clustering procedure 
where each homogeneous region finds another homogeneous region that is closest (in 
measurement space) to it; if the two regions are close enough, they are merged to one 
cluster.' This process of grouping homogeneous regions continues iteratively till there 
are only a few major clusters. 

Determination of Composition of the Clusters 

The composition by terrain categories in the natural clusters produced by 
applying the clustering procedures to the microwave response was performed by a 
manual classification based upon maps and spacecraft imagery. The terrain cate- 
gories chosen were quite gross and mixtures of two single categories (e.g. forest 
or agriculture) were combined to form composite categories (e.g. forest/agriculture 
where there was evidence of both and no predominance of either. 

For some of the clustering exercises, the number of samples was very large 
and manually classifying each observation would have been prohibitively time- 
consuming. For such cases a random sample of up to 40 samples was manually classified 
for clusters which contained over 50 samples. For smaii-size samples, however, 
every point was classified manually. 
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•The manual bssignmenf of targets to terrain categories can often be in error 
because of two reasons: the maps were produced many years ago (5-20 years) and 
the land-use could have been somewhat altered; and examination of a gross- resolution 
map like, for example, the potential natural vegetation map calls for a subjective 
Judgement on the part of the manual interpreter. Spacecraft imagery was used where- 
ever available, but for the cross-track modes of operation only the inner angles in 
the cross-track scan were supported by imagery. Besides, large portions of the space- 
craft Imagery were rendered totally useless because of intervening clouds. 

Results of Clustering 

When considering only the land use and physiographic classifications as de- 
scriptors, the results of the clustering analysis were not encouraging. It would appear 
that factors other than those suggested by major biomes alone should be considered. 

The combinations of sensor responses that were considered in this study may not have 
been optimal for partitioning the terrain according to the physiographic features. Due 
to lack of time only a few clustering exercises were conducted. Future studies should 
consider clustering procedures with variables not considered here and should fry to 
determine the casual relationship between some key parameters of the terrain and the 
microwave response. For example,- one could introduce soil moisture variation as 
a terrain characteristic and assess if the clustered microwave measurements are keyed 
by the moisture variations of the surface. Such an analysis is more complicated than it 
would appear because soil moisture variations can often be masked by vegetation 
canopies. Clustering on some set of terrain descriptors and on some set of microwave 
response should provide similar clustered groups. This choice of terrain descriptors and 
set of microwave responses will be a culmination of many clustering attempts with 
various terrain descriptors and various sensor configuration responses. 

The results reported here represent a first step in the exploratory process and 
suggest that more sophisitication in the definition of the terrain categories is required. 
These clustering analyses were conducted by considering the following combination of 
sensor configuration responses. 
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Figure B~2a. Cross“track contiguous data takes over U .S.A. during SL2 arid SL3 
- missions considered In designing data catalogue. 





IN-TRACK CONTIGUCUS 



Figure B-2b, In-frack contiguous data segments over U.S.A'. during SL2 and SL3 
missions considered In design of data catalogue. 






•1 . The CTC pitch“offset-29® radiometer/scaHerometer“mode dafa including 
both the backscatter and radiometric temperature for vertical polarization. 
Only North American targets are included, 

2, The CTC pitch- offset- 2 9*^ backscatter-only data for vertical polarizationo 
North America, South America and Ocean data are all included o 

3o The ITC (VV) backscattering coefficient for the middle three angles - 
(17^, 33*^, 41°) for all targets which were viewed by all three angles. 

This included some ocean targets and some targets that were not in the 

U.S.A. 

4, The ITC (VV) backscattering data for the middle three angles (17°, 33°, 

41°) for the long Texas~to~Maine pass on Day 253. This pass was also 
clustered by a spatial clustering procedure for comparison. 

Figures B~2a and B-2b show passes from which the ITC and CTC data over North 
and South America were obtained.' 

Olher clustering exercises were also conducted, using ITNC data over land 
and using ITC data for other combinations of sensor configurations than those listed 
above. These showed similar or worse results, although the number of samples to 
be clustered was in some cases not sufficient for any proper inference o These cases 
are not reported. 

Clustering on Backscattering Coefficient and Radiometric 
Temperature at 33° 

The CTC (VV) rad/scat pitch offset mode was exercised extensively over North 
America. Many target points were obtained from oceanic surfaces as well, but these 
were excluded from the data base. K-means clustering using the radiometer response 
and the backscatter response as the two-dimensional clustering variable produced five 
significant clusters. The total number of footprints was greater than 3000 so only a few 
random samples were manually classified for clusters that contained over 100 points. 

A summary of the results is shown in Table B-1 . The categories are described by 
alphanumeric codes; an explanation for these are given atop Figures B-3, B-4 and B-5, 
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Table B-1 


D “ Desert 

R1 - Range - Grassland 
R2“ Range - Thornbush 
R3 » Range - Savanna 


IVU p Range / Agriculrure 
A « Agriculture 
1V12 “ Agriculture /.Forest 
IV\3 “ Range / Forest 


F “ Forest 

T =? Tropical Rain Forest 
W “ Water 
U “ Urban, 

n “ Non VegetatedJerrain 



NUMBER 'Of fOOTPRINTS 

TOTAL 


D 

R1 

R2 

R3 

Ml 

A 

M2 

M3 

B 

Q 

B 

u 

m 

1 

a 

B 

B 

0 

m 

3 

11 

10' 

B 

B 

B 

B 

B 

45' 

2 

B 

3 

m 


B 

1 

23 

B 

13 

B 

B 

B 

g 

4r 

3 

B 

B 

fl 

B 

2 

7 

B 

B 

• r* 

P 

B 


B 

B 

40"' .. 

4 

‘2 

4 

B 

B 

17 

1 

9 


7- 

K 


B 

B 


5 

10 

5 

B 

B 

3 

3 

4 

__ 

15' 

— 


B 

B 


TOTAL 

12 

12 

■ 

B 

22 

■ 

73 

D 

61 

■ 

B 

B 

“ 

^B 


Results Of Clustering Analysis Considering Rad / Scat 'Responses 
From CTC Pitch 29® Operation. 

( * Only A Few Random Samples Were Taken From .Each Cluster . 
The Actual Number Of Points in Each Cluster Was Far Greater. ) 
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Figure B-3 


100 Cluster 



Categories 


'■ % of Category 


D => DESERT 
R1 « RANGE - Grassland 
R2 “ RANGE - Thornbush 
R3 RANGE “ Savanna 
Ml = RANGE/ AGRICULTURE 
A =» AGRICULTURE 
M2 « AGRICULTURE /FOREST 
M3« RANGE /FOREST 
F » FOREST 

T « TROPICAL RAIN FOREST 
W « WATER BODIES 
.U « URBAN 

N « NON VEGETATED TERRAIN 


PROPORTIONS OF MANUALLY CLASSIFIED PHYSIO- 
GRAPHIC AND LAND/USE CATEGORIES IN EACH CLUSTER. 

n VARIABLE? AND EMISSIVITY AT 30^ 

DATA BASE: am rrc r/s DATA 

TOTAL NUMBER OF CLUSTERS: ^ ^ 


VERTICAL INCIDENCE 



Percent' 


Figure B-4 


(O 


D = DESERT Ml = RANGE / AGRICULTURE F == FOREST 

R1 - RANGE - Grassland A = AGRICULTURE T = TROPICAL RAIN FOREST 

R2 « RANGE -Thornbush M2 = AGRICULTURE / FOREST W= WATERBODIES 

R3 » RANGE - Savanna M3 « RANGE / FOREST U « URBAN 


% of Cluster % of Category N » NON VEGETATED TERRAIN 



Categories Categories 


PROPORTIONS OF MANUALLY CLASSIFIED PHYSIO- 


GRAPHIC AND LAND/USE CATEGORIES IN EACH CLUSTER. 

CLUSTERING VARIABLE: a° AND EMISSIVITY AT 30° 


DATA BASF* 

TOTAL NUMBER OF CLUSTERS: 


ALL GTC R/S DATA 
5 


/ 


VERTICAL INCIDENCE 



Percent 



PROPORTIONS OF MANUALLY CLASSIFIED PHYSIO- 
GRAPHIC AND lAND/USE CATEGORIES IN EACH CLUSTEPn. 

CLUSTERING VARIABLE; ^ AND EMISSIVITY AT 30°, VERTICAL INCIDENCE 

DATABASE; ALL CTC R/S DATA . 

TOTAL NUMBER OF CLUSTERS; 5 




Category 

D ~ Desert AF ~ Agriculture - Forest 

R - Range / Grassland RF : Range / Forest 

RA - Range / Agriculture F - Forest 

A - Agriculture 

Cluster Analysis of Cross-Track-Contiguous Microwave Data. 
Clusters Based on Pitch 29° Radiometer and Scatterometer Data. 
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Percent of Category 


Except for cluster 1 -which only had 45 samples all the other clusters had samples 
In excess of a hundred; only 40 samples, randomly chosen were manually classified 
for these clusters. Figure B~3, B-4 and B-5 show a bar graph presentation of the 
proportion of each terrain category in each- cluster. The bar graphs show the ratio 
of a particular category to the total number of samples in a cluster and to the total 
number in that category across all clusters. The first ratio is an indication of the 
composition of the cluster, the second ratio shows a measure of separability. 

These results have been presented another way in Figure B~6, where the 
three-dimensional presentation makes it somewhat easier to see trends. Note that the 
horizontal axis Is arranged approximately in the order of increasing vegetation amount 
and size. Clearly the less- vegetated areas are concentrated in the upper left hand 
of the figure, whereas the forest groups other than agriculture-forest are mostly in 
clusters 1 and 2 . Agricultural terrain is more uniformly spread among the clusters. 

These results are more encouraging for ihe use of microwave sensors than is 
at first apparent, but clearly better categorization of the "ground truth" is called 
for. For instance, the areas labelled desert undoubtedly include some range land, 
and those labelled range would in some cases be called desert by a different categorizer. 
The category containing both range and agriculture includes areas in which the agri- 
culture is dry- land farming that should look very much like range, or even desert, to 
the microwave sensors. The genera! category "agriculture" includes everything from 
dry pastures and freshly plowed ground to densely vegetated stands of corn and milo. 

The former are similar in general dryness and roughness to the rangeland and some of 
the desert, whereas the latter are expected to scatter microwave energy as uniformly as 
a dense forest. Although It is somewhat surprising that "forest" should appear to a 
limited extent in category 5 where desert is so strong, it must be remembered that some 
land is classified as forest when it contains rather sparse stands of scrub trees containing 
little moisture. 

In view of these comments, it seems that more appropriate categories for use 
in ihis study would have been based upon density of vegetation (perhaps the biomass) 
rather than upon more traditional land-use major categories. After all, a dry-land 
bean farmer in central New Mexico has more in' common with his cattle- grazing neighbor 
than with an Iowa corn-and-hog operator, and perhaps the categories used should reflect 
this. Unfortunately, this kind of distinction is more difficult to make from the types of 
land-use maps available for this study, and even from the space photos, than was the 
categorization used. 
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Clustering On Backscafter Response at 33° 

The backscatter data from CTC (VV) pitch offset 29° rad/scat and scattero~ 
meter*only modes was used as a one~dImensiona! clustering variable. The number 
of points was over 7000 and the targets were in North America, ‘South America 
(Brazil) and in the oceans. The number of significant clusters was 10 and a summary 
of the results is provided in Table B~2. Once again only a few random samples were 
classified from the larger clusters. 

Without the radiometer data being included, of the four clusters in which 
there were water bodies, one contained a sizeable number of land targets,' Cluster 5 
also shown In Figure B-7 in a bar graph representation contained some desert range- 
land, and a combination of rangeland, agricultural land and forest. The water targets 
were in the Gulf of Mexico and the land targets were in Texas, New Mexico and 
Arizona. The mean backscatter value for this cluster was -15dB and a verification 
showed that all of the land and ocean samples in this cluster registered very close to 
this mean value. Figure B-7 also shows a bar graph representation of cluster 1 \ The 
non-vegetated terrain category is comprised of a playa surface that appeared on 
imagery to be totally bare. Once again, there seems to be no particular partitioning 
of the terrain categories based upon the scattering data. 

Figure B-8 shows a three-dimensional plot of ihese data like that of Figure 
B-6, The large number of clusters in this case makes discrimination more difficult than 
it was in Figure B-6, Note, however, that the forest categories F, RF, and T are 
concentrated in clusters 7, 9 and 10; the only other categories making a significant 
contribution to this group of clusters are RA, and AF. Conceivably significant portions 
of the RA and RF categories contained either dense crops like maturing corn or 
significant number of trees. Regrettably, it was not possible to check these conjectures 
out' point by point because of lack of time. 

The tropical rain forest Is contained completely in clusters 9 and 10, and it 
is the only constituent of cluster 9. Thus, the ability to distinguish this type of dense 
forest seems assured. Most of the water was clearly distinguished into clusters 4, 6, 
and 8, with rhe exception of the small contribution in cluster 5, The reason for this 
is not known. Also somewhat surprising In view of the results of the radscat evaluation 
is that "Agriculture" appears only in the single cluster, 1 . The only explanation seems 
to be that more of the cells classified as agriculture-only were of the same type in this 
set of data than in the radscat data. 
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Table B**2 


D “ Desert 

R1 “ Range - Grassland 
R2” Range -Thornbush 
R3” Range - Savanna 


Ml p Range/ Agriculrure 
A » Agriculture 
M2 “ Agriculture /Forest 
M3 “ Range / Forest 


F “ Forest 

T ■» Tropical Rain Forest 
W “ Water 
U “ Urban. 

n “ Non Vegetated Terrain 


IS 

NUMBER OF FOOTPRINTS 

TOTAL 


D 

R1 

R2 

R3 

Ml 

A 

M2 

M3 



W 

U 

N 

1 

4 

10 

11 

2 

2 

4 

5 

1 

B 

B 

B 

B 

B 

39 “ 

2 

9 

19 

4 

B 

B 

B 

1 

3 

4 


g 

B 

B 

40 

3 

lii 

m’ 

A 

L 

F 

U 

N 

c 

B 

fl 

0 

N 

B 


4 


B 


B 

B 

g 

g 

g 

B 


40 

B 

B 

40" 

5 

6 

6 

B 

1 

1 

B 

2 

B 

g 

B 

23 

B 

B 

39=- 

6 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

15 

B 

B 

15" 

a 

1 

6 

5 

1 

2 

B 

4 

14 

6 

g 

g 

B 

1 

40* 

8 

i 

B 

B 

g 

B 

B 

B 

B 

B 

B 

40 

B 

B 

40" 

9 


B 

B 

g; 

B 

B 

B 

B 

B 

40 

B 

B 

B 

40" 

10 

g 


5 

Bl 

B 

B 

B 

2 

B 

27 

B 


B 

40" 

TOTAL 

20 

47 

25 

4 

5 

4 

12 

20 

10 

67 

118 

g 

1 

333 


Results Of Clustering Analysis Considering Scatterometer Responses 
From CTC Pitch 29° Operation. 

* 

( * Only A Few Random Samples Were Taken From Each Cluster . 

The Actual Number Of Points In Each Cluster Was Far Greater. ) 
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Figure B-7 


D « DESERT Ml = RANGE / AGRICULTURE F - FOREST 

. R1 = RANGE - Grassland A = AGRICULTURE T = TROPICAL RAIN FOREST 

R2 = RANGE- Thornbush M2 = AGRICULTURE / FOREST W = WATERBODIES 

R3 = RANGE - Savanna M3 = RANGE /FOREST U = URBAN 



D R1 R2 R3 Ml A M2 M3 F T W U N D R1 R2 R3 Ml A M2 M3 F T W U N 

Categories . Categories 


PROPORTIONS OF MANUALLY CLASSIFIED PHYSIO- 

IlXtong categories in each cluster. _ 

DATA BAS^f ALL^CTC^SCA^'^^Ti''^ INCIDENCE 

TOTAL NUMBER OF CLUSTERS: - lo ^ 




Category 


D “ Desert 

R - Range /Grassland (North America) 
RT- Range /.Thornbush (Brazil) 

RS - Range / Savanna (Brazill) ■ 

RA - Range/ Agriculture 
A - Agriculture 


AF - Agriculture/ Forest 
RF - Range / Forest 
F - Forest . 

T - Tropica! Rain Forest 
W- Water 


Cluster Analysis of Cross-Track-Contiguous Pitch 29° 
Scatterometer Data. 
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Percent of Category 


One conclusion that might be drawn from this data set is that the number of 
clusters should be deliberately reduced to something less than 10 by increasing the 
permitted data~space radius. Clearly the sensor is not capable of making real dis“ 
tinctions of this many categories of the type obtained from maps, and allowing the 
number of clusters to remain large tends to confuse interpretation. 

The ITC (VV) data were sorted so that the bpckscatter responses of a target 
to the three angles considered were grouped. This was the three-dimensional measure- 
ment vector used for clustering, A K-means clustering procedure produced 10 significant 
clusters. Three others were discarded because of obvious errors in the computation 
of backscatter or the ephemerides of the target or both. The proportions of the terrain 
categories in each cluster are shown in Table 3. 

One can see that water is always separable from land (clusters 3, 9 and 10), 

The three clusters of water correspond to distinctly different backscattering coefficients 
for each. The probable cause of the three separate water clusters is wind speed. 

Below a wind speed ot 4-6 knots the backscattering coefficient at these oblique incidence 
angles falls very rapidly. 

Figure B-9 illustrates in three dimensions the results of this clustering. The 
clear separation of water is apparent, but it is also apparent that the forested categories 
are concentrated in clusters 2 and 8, whereas clusters 5, 6, and 7 contain mostly agri- 
culture and range. Cluster 4 contains range alone. The amount of agriculture appear- 
ing in cluster 2 is explainable as before is land that contains dense crops, but the pre- 
sence of range in the same category as the forests cannot be explained. 

This data set shows some promise for classification using the scatterometer if 
better categories were developed, but it is not highly encouraging. Perhaps a com- 
bination of better categorization with finer resolution that permits the cells to be more 
homogeneous would work better, but of course this is conjecture. 
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Table B-3 


D “ Desert 

R1 «• Range - Grassland 
R2 " Range “Thornbush 
R3 *> Range - Savanna 


Ml “ Range /Agriculture 
A * Agriculture 
M2 - Agriculture /Forest 
M3 » Range / Forest 


F “ Forest ■ 

T “ Tropical Rain Forest 
W - Water 
U - Urban 

M ° Non Vegetated Terrain 



NUMBER OF FOOTPRINTS 

TfTT ! 


D 

R1 

iR2 

R3 

Ml 

A 

M2 


F 


IV 

U 

M 

1 U 1 r^L 

1 

B 

4 

■ 

■ 

B 

1 

B 

B 

B 


B 

B 

B 

5 

2 

B 

7 

B 

B 

B 

4 

25 


15 


B 

B 

B 

51 

3 

B 

B 

B 

B 

B 

B 

B 

B 

B 


3 


B 

3 

4 

un 

12 

B 


B 

B 

B 

B 

B 


B 

B 

B 

12 

5 

B 

3 

B 


fl 

3 

9 

B 

B 


B 

B 

B 

15 


B 


B 


B 

8 

B 


B 

B 

B 

B 

B 

21 

7 

B 

10 

B 



4 

B 


B 

B 

B 

B 

B 

14 

8 


6 

B 

i 

Bi 

1 

19 

B 

16 

B 


B 


42 

9 

B 

B 

B 


Bi 

B 

B 

B 

B 

B 

• 9 

B 

B 

9 

■M 

B 

B 

B 

i 

B 

B 

B 

B 

B 

B 

15 

B 

B 

15 

TOTAL 

■ 

55 

■ 

■ 

B 

21 

53 

■ 

31 


27 

B 

B 

187 


Proportions Of Land/Use And' Physiographic Categories Clusters 
Produced By Considering Backscatter Coefficients At 42^, 33^ And 
17^ From ITC (VV) Data. Only U. S. And Some Ocean Targets Are 
Considered. 
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Figure B-9 . 
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Category 


R - Range / Grassland 
A - Agriculture 
AF - Agriculture / Forest 
F - Forest 
W - Water 


Ciuster Analysis of In-Track-Contiguous Microwave Data. 
Clusters Based on Multi-Angie Scatterometer Data. 
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Percent of Categor7 


• Comparison of K-Means and Spatial Clustering Using 
Backscattering Coefficienf at 42°, 33° and 17° 

An ITC (VV) pass on September 10, 1973 which viewed regions from Texas 
to Maine (see Figure B~l) was considered as a sample where two clustering approaches 
were compared. The middle three angles (42°, 33° and 17°) backscatter response 
was considered as the three-dimensional variable, K-means clustering results were 
obtained from the clustering exercise described above. Only clusters with a signi- 
ficant number of samples were chosen. The clusters which contained only one or 
two samples in a particular category and for which the targets were from this sample 
pass were ignored because they would be misleading if shown. The three clusters 
produced are shown in Table B-4.' 

The backscattering coefficient at each of the angles considered was then quan- 
tized separately by an equal probability quantizing algorithm. This three angle back- 
scatter was then considered as the three-dimensional variable for spatial clustering. 

Spatial clustering has been used primarily on images, therefore to make our data appear 
in an image- like format we repeated the variables n times to create n columns of the 
image . 

Two approaches with spatial clustering were attempted: for one, the gradient 
threshold was set at unity and regions of homogeneity were first found before the 
clustering procedure was applied; and two, each row of the psuedo- image (corresponding 
to each original measurement vector) was considered to be an independent quantity 
for clustering. Note that in the first case regions of transition (when the gradient thresh- 
old has been exceeded) would not register and those data points would not be considered 
in the subsequent clustering procedure. 

Tables B-4, B-5 and B-6 show the results for the K-means clustering, spatial cluster 
ing with each element considered In the clustering and with a gradient threshold set 
at unity respectively. Although the ratio of the elements in the columns for agri- 
culture, agriculture/forest and forest vary there is no clean-cut differentiation of terrain 
types possible with any of the three. 
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D ” Desert Ml " Range /Agriculrure ■ F * Forest 

R1 ” Range - Grassland A ~ Agriculture T " Tropical Rain Forest 

R2- Range - Thornbush M2 - Agriculture /Forest W ■ Water 

R3 - Range - Savanna A0 - Range / Forest ■ U ■ Urban 

» n “ Non Vegetated Terrain 



Table B-4 K-means Cluster Analysis Using o° at 42°, 33° and 

17° From ITC (VV) Mode for Texas-to-Maine Pass on Day 253 



Table B-5 ^®sults of Spatial Clustering Analysis for Same Data 
Set Used for Table B-4. . No Gradient Threshold. 



Table B-6 R^^^ilts of S'patiai Clustering Analysis for Same Data 

Set Used for Tables B-4 and B-5. Gradient Threshold = 1. 

24 ■ 

































































Conclusion 


Based upon the results of a few clustering exercises using certain combinations 
of Skylab radiometer/scatterometer responses, if appears that the terrain features 
identified by maps and imagesy are not necessarily separable in their microwave re~ 
sponse. It may very well be that the microwave response to other combinations of 
sensor configurations could produce results which can allow a separation of terrain 
based upon the microwave response. The Skylab radiometer/scatterometer has a 
gross resolution and the assignment of observations into terrain categories was based 
upon examination of maps (supplemented by spacecraft imagery) which calls for a 
subjective decision. Furthermore, the use of land-use rather than vegetation-density 
categories was probably a mistake. These may have been factors which contributed 
to the results reported here. 

Variation in vegetation biomass and particularly variations in soil moisture 
within and among categories could have been a factor, but these are not obtainable 
from the ground information available. 

Clearly the microwave responses do separate into different clusters. Proper de' 
scriptlons of the content of these clusters would be extremely interesting, but this 
must be deferred until later e>qDeriments for which more extensive ground information 
can be collected specifically for the experiment. These as yet undescribed categories 
may be very useful for some applications, but apparently gross land-use determination 
is not one of the suitable applications. 
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APPENDIX C 

BACKSCAHER RESPONSE AT 13,9 GHZ FOR MAJOR 
TERRAIN TYPES AS SEEN FROM ORBIT 

A. SobH/ R.K. Moore, andS.T, Ulaby 
RSL TR 243-4 

ABSTRACT 

The S-1 93 radar scatf-erometer that flew aboard Skylab provided information on 
radar returns from many kinds of terrain. Many of the "footprints” of the scatferometer 
have been classified into terrain types, primarily in terms of land use, by examination 
of physiographic and land-use maps along with Skylab photographs where clouds did not 
interfere. The mean scattering coefficients for vertical polarization observed at 
angles of incidence between 1 .5° and 46*^ are reported here for water bodies. The 
standard deviations of the measurements are also reported. Means vary from about 13 dB 
at 1 ,5° for the salt flats of Utah to about -18 dB for water at 46°, For more normal 
land categories they range from about +2 dB at 1 .5° for desert to -12.5 dB for rangeland 
at 46°. Since the categories are broad and not based upon expected microwave 
response, but rather on physiographic and land use descriptions, the ranges between 
upper and lower standard deviations overlap for most categories. 

INTRODUCTION 

The space operation of the S-1 93 scatterometer, operating at 13.9 GHz 
(2.16 cm wavelength) on board Skylab provided data over many regions of the world. 
Statistics of the backscatter response from various regions of the world have been 
reported elsewhere. Comparisons are made here of the backscatfering coefficients 
measured for various terrain categories which were identified by an examination of 
physiographic and land use maps and spacecraft imagery (from S190A and S190B 
cameras on board Skylab). The terrain categories chosen are rather gross because the 
target resolution size of the S-1 93 scatterometer is large (ranging from 100 square 
kilometers to 400 square kilometers) . Since this effort involves comparisons of 
backscatter response over "footprints" that were classified by manual interpretation, 
the data base consists of only a small sample from the many passes made by the S-1 93 
scatterometer. The emphasis in selecting candidate passes has been on finding areas 
where the land use and physiographic features were easily accessible and where large 
areas of homogeneous terrain were known to exist. The large extent of rain forests in 
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Brazil is a prime example of the latter case , The backscatter response from ocean 
surfaces, computed by considering many passes over the ocean, is shown for comparison. 
The data presented here are for vertical polarization only. 

THE SKYLAB RADAR SCATTEROMETER EXPERIMENT 

One of the sensors observing the earth from Skylab was a radar scatterometer 
(termed the S-193), operating at 13.9 GHz. Characteristics of this instrument have been 
described in various NASA pubi icatlons and in some Journals, so only the briefest 
summary will be included. 

The Skylab scatterometer was part of a composite radiometer/scatterometer and 
altimeter all operating at 13.9 GHz and sharing the same parabolic antenna and much 
of the rf hardware. The beamwidth (two-way) of the mechanically scanned dual- 
polarized antenna was 1 .45°. The scan of the antenna was programmed and selection 
of one of four possible modes was made by the astronaut. 

1 . In-Track Non-Contiguous (ITNC — overlapping measurements at angles of 
0°, 15°, 29°, 40° and 48° from nadir, with 100 kilometers spacing between each set of 
targets with measurements at all five angles). 

2. Cross-Track Non-Contiguous (CTNC — measurements at the same angles of 
incidence as ITNC but perpendicular to the track so that there Is a 100 kilometer 
spacing between every target point and no overlap). 

3. In-Track Contiguous (ITC — measurements at the same angles as the non- 
contiguous modes but in pitch (along track) and at a spacing of 25 kilometere between 
adjacent targets with all five angles). 

4. Cross-Track Contiguous (CTC — measurements within +1 1°, perpendicular to 
track about a scan initialized point. The scan could be initialized at nadir, or 15°, 

30°, 40° in pitch or roll but not both) . 

The non-contiguous modes were exercised mainly over the ocean although some 
were also flown over land. For this study the contiguous modes were the major sources 
of data over land and the non-contiguous modes were the exclusive sources for data 
over the ocean; The non-contiguous modes provided coverage with all four transmlt- 
receive polarization pairs while the contiguous modes provided coverage with any 
selected polarization pair. The target size varied with incidence angle and ranged from 
a minimum of 97 square kilometers (a circle with a diameter of approximately 11 kms) to 
to a maximum of approximately 387 sq. kms (an ellipse 17.3 by 28 kilometers). 
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Data from the first and second occupancies of Skylab were used in this study. 

The passes chosen over land are shown in Table 1 . These passes are predominantly 
over the continental U.S.A. although there are a few over South America. There were 
a total of 10 passes with the ITNC mode configuration over ocean surfaces which were 
used to compute the statistics of the backscattering coefficients over the ocean. 

The major terrain categories into which all terrain was to be classified were: 

1 . A — Agriculture 

2. D1 — Desert (Nevada, Arizona) 

3. D2 — Desert (New Mexico, Texas) 

4. F — North American Forest 

5. R — Rangeland 

6. SF — Salt Flats, Utah 

7. SH — Sand Hills, Nebraska 

8. TF Tropical Forest, Brazil 

9. W — Water, Ocean 

The comparisons were made by considering a number of samples for each polariza- 
tion and angular group. 

CHARACTERISTICS OF BACKSCATTER FOR VARIOUS TERRAIN 

For terrain categories with sufficient samples in a particular 2° incidence-angle 
group centered at 1 .5°, 17^, 33^ and 46^, histograms of the distribution of the back- 
scattering coefficient were generated. Along with these histograms, standard deviations 
of the mean backscattering coefficient were also computed. Where the sample size for 
a particular terrain category and incidence angle group was insufficient to warrant a 
histogram, only the mean value is provided. Some terrain categories had enough 
samples at each of the four angles considered to produce histograms, while others only 
had sufficient data at one angle. 

The angles chosen for describing the backscatter coefficient are dictated by the 
availability of data. These angles are the true incidence angles for the in-track modes 
of operation; the cross-track contiguous mode with the appropriate pitch- or roll-offset 
also provides data at or near these angles. Over land, the. number of samples ranges 
from a maximum of 99 for Forest at 46° to a minimum of 8 for .Desert (Nevada, Arizona) 
at 1 .5®. The number of samples over the ocean surface always exceeded 82 for all 
four angles of incidence considered. 

The backscattering coefficients for the various terrain categories are presented 
for each of the four incidence angle groups separately. Figure 1 shows the distribution 
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of backscattering coefficient for the different categories for an incidence angle 
centered at 1 .5°. The range between the upper and lower values that are one standard 
deviation from the mean at this angle are greater than for any of the other angles for 
every surface. The range between these upper and lower bounds for the ocean are 
greater at all angles than for any land category. The mean for the ocean is lower than 
that for the Salt Flats in Utah but the upper bound (hereinafter the one sigma value 
is called the bound) for ocean is higher than the mean for the Salt Flats. As an example, 
the histogram of the distribution for the ocean for vertical polarization and for the angle 
group considered here is shown in Figure 2. The mean value for cr° (the backscattering 
coefficient) is 12.32 dB. This is about 10 dB higher than the backscatter over an 
ensemble of targets over North America. The mean value for Utah Salt Flats (13.5 dB) 
is quite expected because the Utah Salt Flats appear like a smooth reflector at 2 cms 
wavelength. 

The backscatter responses for the remaining five terrain categories which 
registered data for this angle (see Figure 1) show a great amount of overlap. Rangeland, 
Forest, and Agriculture had enough samples to produce histograms. The spread between 
the upper and lower bounds for all three of these groups is approximately 7 dB. The 
range between upper and lower bounds for these categories overlap the means for desert 
in Nevada, Arizona, and the sand-hills region in Nebraska. In fact, the lower bound 
of responses over the ocean falls within these ranges. The difference in the means for 
these land categories is very small . 

At 17° incidence (Figure 3) the range between the upper and lower bounds 
becomes smaller for all terrain categories. The means for all terrain categories are 
much lower than at 1 .5° (as expected) . The mean backscatter over the ocean still 
registers higher than land categories. The difference between the means for forest and 
for rangeland at 1 .5° was approximately 1 .9 dB, at 17° the mean for the two categories 
is almost the same. 

The maximum number of data included in this ^udy is for an incidence angle 
group centered at 33°, with seven of the nine categories containing enough samples to 
permit the generation of histograms. There is a reversal in the comparative level of the 
ocean backscatter (Figure 4). The mean for the ocean at 33° is now seen to be roughly 
5 dB lower than most land targets. There is a wide spread between the upper and lower 
bounds (8 dB) for ocean; the corresponding spread for all land targets is" much 
less. Once again the means for all land targets except the tropical forest in Brazil lie 
within or very close to the bounds between upper and lower values for all land categories. 



The tropical rain forest with its dense vegetation canopy seems to appear rougher to the 
scatterometer and registers a backscatter higher than the corresponding North American 
forest . 

Less data exist at 46° than. at 33°, The range between upper and lower bounds 
for the ocean surfaces is seen to be very large {10 dB) . This range is probably larger 
than would be expected for only a 2° Interval of incidence angle around 46°. To 
increase the sample size, incidence angles from 42° to 50° were included in this 
angle group. The mean for the oceanic backscatter is now seen to be some 10 dB below 
that for agricultural terrain and about 7 dB below that for rangeland. The range between 
upper and lower bounds for rangeland is quite large (9 dB) . Surprisingly the mean for the 
agricultural terrain is higher than that for forest and rangeland. Perhaps, the small range 
between upper and lower bounds for agricultural terrain must be due to small sample 
size (48 samples). 
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TABLE 1 . S-193 DATA PASSES CONSIDERED 



Mode/Angle 

Area 

Description of Area 

156 

CTC POO 

Ui^ah 

Salt Flats and Great Salt Lake and 
Desert 

156 

CTC P29 

Texas 

Woodland, Agriculture, Grassland 

161 

CTC P29 

Brazil 

Mostly Tropical Forest 

162 

CTC P29 

Oregon 

Agriculture, Rangeland and Forest 

162 

CTC P29 

Venezuela, Brazil 

Mostly Tropical Forest 

165 

CTC P29 

Brazil 

Mostly Tropical Forest 

165 

CTC P29 

Nevada, Arizona 

Mostly Desert 

215 

CTC P29 

Nevada, Utah 

Forested Mountains, Rangeland and 
Portions of Salt Lake 

220 

CTC P29 

New Mexico, Texas 

Woodland, Agriculture and Rangeland 

221 

CTC POO 

North Dakota 

Agriculture and Range 

249 

CTC POO 

Colorado, Nebraska 

Rangeland, Agriculture 

250 

CTC P29 

S.W. Mexico to New 
Mexico 

Mostly Desert 

253 

CTC POO 

Colorado, Nebraska 

Forest, Agriculture and Sand Hills 

223 

ITC 

Nevada' 

Mostly Desert, Little Forest and 
Mountains 

250 

ITC 

New Mexico to Kansas 

Agriculture and Forest and Rangeland 

253 

ITC 

Texas to Maine 

-Farmland, Rangeland, Forest, Some 
Lakes 

217 

ITNC 

Washington Coast to 
Idaho 

Agriculture and Forest 

217 

ITNC 

Idaho to Oklahoma 

Agriculture, Mountains, Sand Hills, 
Rangeland 
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COMPARISON OF SKYLAB SCATTEROMETER 
DATA WITH PRIOR MEASUREMENTS 

King and Moore^ conducted a survey of terrain scattering measurement 
programs. The categories and angles for which they provided data during the pre- 
Skylab era were selected so that a ready comparison with data obtained from the 
Skylab scatterometer could be made. Due to the extensive effort involved in 
locating and identifying homogeneous targets by imagery and maps, only a few 
categories could be represented by the Skylab data. The incidence angles which 
King and Moore used in their comparisons were the design values of the incidence 
angles expected during the space operation of the Skylab scatterometer. The 
incidence angles actually achieved during the space operation deviated slightly 
from these design values. This should have minimal effect for the angles away from 
nadir but for the angles near vertical, one must allow for a variation due to inci- 
dence angle. The pre-Skylab data have been reported at 0*^ whereas the Skylab 
data are reported at 1 ,5°. 

King and Moore reported on both vertical and horizontal polarization. In 
this study we have considered only vertical polarization. The fact that the response 
for the two polarizations is quite similar and correlated allowed us to spend more time 
getting a larger data base for just one polarization. The categories considered by 
King and Moore are somewhat different from ours, e.g., they make no distinction 
between types of forest whereas we find that Brazilian Tropical Forest formations 
show decidedly different backscatter characteristics than any of the less dense 
North-American forests. The categories used by King and Moore are necessarily 
designed so that the bulk of backscatter data collected by programs (aircraft-based 
and ground-based)' could be reported. For example, categories like road surface, 
which could not possibly show up in the gross resolution of the Skylab scatterometer, 
were reported. 

We have taken the data from King and Moore's report for the categories for 
which we could locate Skylab scattering coefficient data. The presentation of the 
comparisons Is given In a manner similar to theirs, i. e., in a bar graph representing 

if 

the bounds (where available) and the mean values. The bounds as reported by King 
and Moore are also one sigma spreads. 

*"Bounds" Is used as in the preceding section to mean values lying one standard 
deviation from the mean. 
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Figure 6 shows a comparison of the backscatter response at nadir with the 

Skylab scatterometer and measurements by other programs. Naval Research Labora- 

2 

tories conducted two experiment programs, one from aircraft-based sensors and the 

3 

other from bridge-mounted sensors . The frequency for which the comparison is 
shown with Skylab data is X-band (8.9 GHz). NRL did not provide bounds for their 
data for grassland; their average value estimated at nadir is approximately -12 dB. 

This is much lower than the Skylab-obtained value. Sandia Corporation^ conducted 
measurements near vertical at 3.8 GHz. A comparison can be made between Sandia 
data over farmland,, and data obtained from a fan-beam scatterometer operating at 
13.3 GHz on an aircraft over agricultural terrain in Kansas, with the Skylab data. 

It is seen that the bounds all overlap, but, the mean for the Skylab data is a little 
lower, and, the range between upper and lower bounds for the Skylab data is much 
less than the ranges shown for the finer resolution sensors. There is a great variation 
between different types of farmland, but with the great spatial averaging involved 
in the Skylab scatterometer measurements, the' dynamic range seems to have been 
reduced. 

The response over forests is reported by NRL, land and air, and the difference 
between the mean values for these two programs is about 30 dB for the same target 
category. Sandia Corporation’s measurements (3.8 GHz) show a mean value and range 
comparable to that for the Skylab radar. The forests In Sandia's measurements are 
of two densitites (Maine and Minnesota). The Skylab radar viewed regions which 
were registered on the topographic maps as forests but no density of the foliage 
can be estimated. A comparison of the mean desert response at 0° to the NRL air- 
borne measurement shows a difference of 4dB. 

The comparisons at around 16^ are given in Figure 7. Once again the mean 
over grassland obtained by NRL is much lower than that from the Skylab scattero- 
meter. Over farmland, we find again that the dynamic range of values from the 
Skylab radar is much smaller than reported by aircraft end land-based measurement 

5 

programs. There seems to be a variation in the means as well . The Ohio State data 
is at X-band. The resolution cell size for these data is of the order of one foot on a 
side! Notice the extremely narrow range between upper and lower deciles for the 
forest category for the Skylab data. The mean is lower than those measured by NRL 
and Sandia. The desert response from the Skylab radar compares very well with that 
reported by NRL (airborne) measurements. 

Comparisons at 31° are shown in Figure 8. As usual NRL measurements over 
rangeland are over lOdB lowet in the mean than the Skylab data. A significant 
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Figure 6. Comparison of Skylab Results with Previous Aircraft and Ground Scattering 
Measurments. Angle of Incidence 0°. 
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Figure 7. Comparison of Skylab Results with Previous Aircraft and Ground Scatteri 
Measurments. Angle of Incidence 16 . 
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Figure 8. Comparison of Skylab Resuifs with Previous Aircraft and Ground Scattering 
Measurments, Angle of Incidence 31°. 





Figure 9. Comparison of Skyidb Results with Previous Aircraft and Ground Scattering 
Measurments. Angie of Incidence 43*^. 







reducMon in fhe dynamic range can be seen in the backscatter response measure- 
ments over agricultural terrain. The range between upper and lower bounds for 
Skylab data is only 3dB. This compares with ranges in the neighborhood of lOdB 
for the finer-resolution measurement programs. Once again forest has a very small 
range between upper and lower bounds (I dB). Desert has a mean value comparable 
to that of the NRL measurements but a spread which is much more than that for 
forests . 

Figure 9 shows a comparison of 43°; the ranges between upper and, lower 
bounds for rangeland is larger at this angle than at any other angle. The mean is 
again approximately lOdB higher than the NRL value. This implies that if one dis- 
counts the 10-or-so dB bias between the two measurement programs, the trends of 
the backscatter response are quite similar for the NRL and Skylab data. The range be- 
tween upper and lower bounds for agriculture Is only 5 dB; this compares with ranges 
in excess of lOdB for the finer-resolution sensors. Forest at 43° had twice as many 
samples of Skylab data as for the other angles; its range between upper and lower 
bounds is 5dB. The mean falls between the NRL land-based and aircraft-based values 
for forest, which show a difference of over 20 dB. Desert had a mean comparable to 
the NRL values . 

SUMMARY OF SKYLAB SCATTEROMETER RESPONSE 
FOR TERRAIN CATEGORIES 

For the terrain categories considered and for almost all angles of incidence, 
the range between upper and lower one sigma values for the Skylab scatterometer data 
is less than the corresponding range measured by the finer-resolution aircraft-and 
land-based sensors. There is a wide discrepancy between the mean values among the 
aircraft and land-based sensors for each terrain category at each angle. The Skylab 
data usually fall within the upper and lower values presented by the prior measurements. 
For rangeland and desert, where only one prior measurement program has reported data, 
the Skylab data show a difference of 10 dB for the former and are comparable for the 
latter. 

There Is a great deal of overlap in the backscatter response from the various 
categories. Figure 10 shows the bounds between upper and lower deciles and the mean 
values for various terrain categories for angles of Incidence from 0° to 47°. The cate- 
gories included in this composite figure are those for which enough samples existed for 
statistical validity. It can be seen from the figure that the range between upper and 
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lower deciles for the ocean is much wider than all the land categories. There is a 
difference between some land and the ocean targets below 10° incidence and again 
beyond 30°. The region between 10° and 30° seems to have a total overlap between 
the response from all kinds of terrain. This would imply that a radar designed for 
discriminating/identifying terrain would not be very capable between 10° and 30° 
of incidence. The response below 10° incidence shows that, apart from regions like 
the Utah Salt Flats which appear like mirror reflectors at 2 cms wavelength, land 
targets register a backscatter response lower than the ocean surface. The tropical 
rain forests in Brazil show the highest backscatter of any targets at 33°. This is in 
accordance with the expected response from rough surfaces. Notice that the range 
between upper and lower deciles for the tropical forests in Brazil and the forests in 
North American do not ever overlap. Surprisingly, the response from desert surfaces 
is quite high at the greater angles of incidence. 
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Figure 10. Comparison of angular backscatter response from various categories. Ranges, where shown, 
are between upper and lower deciles, not standard deviation. 
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ABSTRACT- 

'Microwave data acquired over the Great Salt Lake Desert area by sensors aboard 
Skylab and Nimbus 5 indicate that the microwave emission and backscatter were -strongly 
influenced'by contributions from subsurface layers of sediment saturated with brine. This 
■phenomenon was observed by Skylab's S-194 radiometer operating at 1 .4 GHz, S"’193 
' "'RADSCAT (Radiometer~Scatterometer) operating at 13.9 GHz and the Nimbus 5 ESMR 
(Electrically Scanning Microwave Radiometer) operating at 19.35 GHz, The availability 
of ESMR data over an 18 month period allowed an investigation of temporal variations. 

’ Aircraft 1 ,4 GHz radiometer- data acquired two days after one of the Skylab passes 

confirm the satellites' observations. ESMR data reveal similar responses over the Bolivian 
deserts, which have geologic features similar to those of the Utah desert. 




1.0 INTRODUCTION 


The microwave emissivity, e , and backscafiering coefficleni, of terrain 
surfaces ere functions of the dielectric properties and surface roughness (relative to 
the wavelength) of the ground. The dielectric properties are in turn strongly influenced 
by the soil moisture content. Microwave observations of soil surfaces by active [Ulaby, 
1974a, b; Uloby et d!., 1975a] and passive [Schmugge et al 1974; Newton et al., 

1974; Eagieman and Ulaby, 1974] sensors indicate a high degree of sensitivity to soil 
moisture variations. Due to the nature of the scattering and emission phenomena, the 
scattering coefficient exhibits a positive correlation with soil moisture content, whereas 
the emissivity (and hence brightness temperature) decreases with soil moisture contenf;*- 
In both cases lUlaby et al ., 1975a; Newton et a! ., 1974] longer wavelengths have 
been observed to yield more satisfactory results (in terms of sensitivity to moisture 
variations) simply because, for a given terrain surface, the effects of surface rough” 

-ness on the microwave response ^ackscatter and emission) are reduced as the wave” 
length is made longer since the surface would appear electromagnetically smoother. 

Brightness temperature data acquired by Skylab and Nimbus 5 microwave 
-•^radiometers over Utah indicate a consistent difference in temperature between the 
-Great Salt" Lake Desert area and neighboring land surfaces. The Skylab microwave 
• sensors include a 13.9 GHz Radiometer”Scatterometer (RADSCAT) designated S-193 
and a L”band radiometer operating at 1 .4 GHz designated S~194. The microwave 
sensor aboard Nimbus 5 is a 19.35 GHz Electrically Scanning Microwave Radiometer 
(ESMR). 

During Skylab Pass 5 on June 5, 1973, S“193 measured brightness temperature 
values as low as 200°K for some parts of the Great Salt Lake Desert in comparison to a 
270°-280°K range for areas outside the desert. In conjunction with the low temperature 
values, the measured scattering coefficient of the same general area was more than 15 dB 
higher than the scattering coefficient of areas outside the desert. Similar observations 
to those indicated by the S-193 radiometer were also evident in the data acquired by 
S-194 and ESMR for the same date. Another Skylab pass on August 8, 1973 and 
numerous ESMR passes over an 18 month period confirm that these observations are 
in response to specific characteristics of the Great Salt Lake Desert soil material. 
Detailed analysis of the data and the hydrology of the region has led us to believe that 
a significant contribution to the measured emitted and backscattered energy is from 
subsurface layers of brine. 



As will be discussed in fhis paper, fhe availabilify of microwave data over 

the Great Salt Lake Desert at three different wavelengths has proven very useful in 

• 

the analysis and interpretation. Moreover, the scatterometer data at 13.9 GHz has 
also served to complement the observations made with the passive sensors. The daily 
global coverage of the Nimbus 5 ESMR make possible observation of other salt deserts. 
In particular, low brightness temperatures have been observed over the salt deserts of 
the Aiti Plano region of Bolivia. 

2.0 DEVELOPMENT AND CHARACTERISTICS OF THE GREAT SALT LAKE DESERT 

The test site under investigation trends in a northwest-southeast direction across 
the Great Sait Lake Desert (Figure 1). The narrow side represents the coverage on the 
• ground by S-193 RADSCAT as the antenna was scanned in the cross-track-contiguous 
■mode at 0° forward pitch during.'The Skylab June 5, 1973 descending pass (northwest to 
.^southeast direction) over Utah. S-194 coverage comprised approximately 74 per cent of 
the test site and complete coverage by ESMR is available for numerous passes. Two 
lines are shown in Figure 1, a horizontal line representing a transect across the desert 
ot 40° 30' North latitude, the elevation profile of which is shown in Figure 2, and a 
. NW-SE line indicating the coverage by the'NASA airborne 1.4 GHz radiometer flown 
on August 10, 1973. Several types of terrain are also shown including lakes (Great 
Salt Lake and Utah Lake), saitflats, mud flats, mountains, and the Great Salt Lake 
Desert. 

2.1 Hydrologic Properties 

fifteen to twenty thousand years ago, coincident with the beginning of the 

■retreat of the most recent glacial advance over the northern portion of the North 

American continent, a large area of northwestern Utah was covered by ancient Lake 

Bonneville which, through a complex history of contractions and expansions, ultimate- 

2 . ' ■ 

ly was reduced to the 4270 km now occupied by the Great Salt Lake. To the 
southwest of the lake and connected to it by a narrow threshold is the Great Salt Lake 
Desert (Figure 1) , approximately 9 m above the present level of the lake. Until 
approximately 10,000 years ago [Eardley, 1962], this area was covered by the waters 
of Lake Bonneville into which were washed the sediments now composing the lake bed 
"deposits which underlie the salt encrusted surface . These deposits underlie the major 

portion of the Great Salt Lake Desert bglow altitudes of approximately 1300 meters 
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ELEVATION ABOVE MEAN SEA LEVEL IN METERS 


ELEVATION PROFILE ALONG A TRANSECT AT 40° 30' NORTH LATITUDE 



Figure 2, 


Elevation profile along a transect at 40°30* North latitude. 



and are predommai’ely clay and sill’ wilh variable sail content' [Stephens, 1974]. 

From this bleak terrain surface rimmed by mountains (Figure 2) several areas show 
deviation: to the west in the general vicinity of Wendover, Utah lay the Bonneville 
Saif Flats and along the east side of the desert are gypsum sand dunes which have 
developed through the ablation of the desert surface by the westerly winds. In 
addition, isolated mosses of bedrock protrude upward from the central and marginal 
areas of the Great Salt Lake Desert. 

Lake bed clays and silts and crystalline salt form a shallow brine aquifer a 
maximum of 7.6 meters thick which covers the bulk of the surface of the Great 
Salt Lake Desert (Figure 3), only the northern portion of which has been studied in 
detail [Stephens, 1974], 

Generally, in the area where lake clays form the surface, the depth to brine 
is estimated to range between 60 and 90 cm near the center of the desert floor and 

2.1 “ 2.7 m at the margins [Nolan, 1928], although capillary action in the fine 
grained sediments may raise the water in excess of one”ha!f meter above the water 

i-table which causes the surface to remain perpetually moist. In the Bonneville Salt Fiats 
-although the surface ot the salt bed is rigid, the salt remains saturated v/ith brine 
to within a few inches of the surface [Stephens, 1974], the water table ranging 
--between approximately 15 and 20 cm below the surface. With the precipitation in 
-the central part of the desert averaging less than 13 cm per year, the possibility of 
.^accumulating standing water on the desert surface or of elevating the water table ' 
-above its normal 60"90 cm position below the .desert's surface is minimal . Runoff 
from the highest parts of 'the mountairs in the area averages less than 2.5 cm 
-annually [Badgley et al., 1964: Busby, 1966] and runoff during the brief periods of 
—rapid snow melt generally infiltrates the stream channels downslope and only locally 
•spreads out over the desert floor, most frequently in the Bonneville Salt Flats area. 
Recharge is by infiltration of precipitation and lateral subsurface inflow, the brine 
moving through the lake beds by intergrannular flow through layers of salt Impregnated 
clay and through open joints (Figure 4). 

2.2 Dielectric Properties 

Except for d few isolated segments, the Great Sait Lake Desert is characterized 
by a very flat surface. Hence, as a first order approximation, the specular surface 
-model can be applied to calculate the emissivity of the' surface in terms of the power 
reflection coefficient R [Moore, 1975]: 
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NORTHERN PORTION OF THE 
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Figure 3, Northern portion of the Great Salt Lake Desert 
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INFERRED SUBSURFACE STRATIGRAPHIC RELATIONSHIPS NEAR WENDOVER 



Figure 4, Inferred subsurface sfrafigraphic relaHonships near Wendover, Ufah. 



e = 1 - R 


( 1 ) 


Af- nadir, R f-akes ^he form: 

R = 

where k is fhe relaHve complex dielectric constant of the ground: 

• k = k^-jk2 (3) 

At microwave frequencies, typical values of k^ for dry soil lie in the range 
2.5 <k-| ^ 3.5 and ^2 ^ ^ [Hoekstra and Delaney, 1974] leading to an 
emissivity range at nadir of 0.95 < e < 0.91 . In contrast, and k 2 of water 
are much larger in magnitude and are, in general, functions of frequency, temperature' 
end salinity [Stogryn, 1971], Shown in Table 1 are calculated k and e .values for 
pure water, saline water with a salinity of 150 °/oo*and saline soil typical of the 
Great Salt Lake Desert brine layer, all at 23°C . According to Stephens [1974] the 
salinity of the water under the desert floor is generally 150 °/oo or more and the 
--.overage porosity of the soil is 45 per cent. Thus the brine layer is a mixture consist- 
ing of 45 per cent saline water and 55 per cent soil . Due to the lack of dielectric 
.Inconstant data for such a mixture, the entry in Table 1 for the dielectric constant of 
saline soil was calculated using the dielectric values calculated for dry soil and saline 
-water [Stogryn, 1971], each weighted by its respective proportion of the combined mixture. 

Most of the brightness temperature data presented in the next section was 
• oequired on June 5, 1973. The ground temperature at the time of the Skylab and 
-.Nimbus passes was estimated to be 296°K {23®C) [NASA, 1974] , hence ail values 
•shown in Table 1 were calculated at 23°C. 

Whereas the effect of salinity on the emissivity of water is negligible at' 

13.9 GHz and 19.35 GHz, Table 1 shows that the emissivity of saline water at 
1,4 GHz is less than half the emissivity of pure water. Moreover, the calculated 
emissivity of saline soil (as defined above) is smaller than the emissivity of pure waterl 



* 



stands for "ports per thousand". 
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Table 1 . Calculafec) dielec^rIc conslanl" t< and omtsstvUyeof dry soil, pure waher, 
saline water and saline soil at 23°C, 

Dry Soil , Pure Water Saline Water* Saline Soil** 


Frequency 


^d 

k 

P 

"P 

k ^ 

. s 

k 

ss 

'^ss 

1.4 GHz 

3-jO 

0.928 

78.63-J5.57 

0.364 

41.9-J237.4, 0.179 

20. 5“j 106.8 

0.256 

13.9 GHz 

3~j0’ 

. 

0.928 

52.4-j35,7 

0.382 

31.4-J40.6 0.394 

15.7-J18.3 

0.524 

19.35 GHz 

3-jO 

0.928 

40.4-J37.2 • 

0.396 

25.7-J35.6 0.416 

13.2-J16.0 

0.549 


*Salinity = 150 °/oo 

**k =0.55 k. +0.45 k 
ss d s 











3.0 MICROWAVE OBSERVATIONS 


Al- sateMUe alHtudes, the measured brightness temperature of a flat surface 
Is given by [Moore, 1975]: 


*'’a ^ ^BS ^ ^ ^d^ ^u 


(4X 


where is atmospheric transmittance, Tg^ is the brightness temperature of the 
surface, is the brightness temperature of the total downward radiation (atmospheric 
plus cosmic) and is the brightness temperature of the upward atmospheric radiation. 
Igs Is the product of e and the ground thermometric temperature T^. Under clear sky 
conditions, atmospheric effects are negligible for the dry desert atmosphere of the 
Great Salt Lake Desert area, particularly at 1 .4 GHz and 13.9 GHz, After calculating 
T , and Tj for the June 5, 1973 atmospheric conditions, the error between the 
measured brightness temperature Tg and the brightness temperature of the surface ^BS* 




( 5 ) 


•was estimated by calculating e using Eq. 4 for the range of values of Tg measured by 
the satellite radiometers. The following results were obtained: 


! ^2,3°K 
<2.5°K 
<9,1®K 


at 1 .4 GHz. 
at 13.9 GHz 
at 19.35 GHz 


Similar results (within 0,2°K) were also obtained for the August 8, 1973 Skylab and 
Nimbus 5 passes over the Great Salt Lake Desert, In addition to these two passes, data 
reported herein includes Nimbus coverage over a period of 18 months during which the 
data acquired on some passes were more severely influenced by atmospheric conditions. 
During the drier winter months, calculations for the mid“!atltude winter standard 
atmosphere profiles yield a smaller difference, AT 5.9°K. The value of the dif- 
ference for the Bolivian salt desert would be less than 5°K due to the higher altitude, 
and reduced water vapor content. 

If we model the Great Salt Lake Desert as a stratified medium consisting of a 
relatively dry surface layer of height h (layer 1) covering a substrate consisting of 
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saline soil (layer 2) with bofh layers having ^he same I'hermomefric ^empera^ure T 
the brightness temperature of -the surface can be expressed as (adapted from King 

2 


?1970]): 




) ( 1 - 


( 6 ) 


where R is the power reflection coefficient at the interface between two semi~ 
infinite homogeneous media. The subscript "a" refers to the air medium above 
laye'" 1. The effect of the attenuation by layer 1 is accounted for by the 
transmittance T-jt- 

■ T^(h) = exp(-2oh) (7) 


where 2a is the power attenuation coefficient of layer 1 . 

In the absence of the top layer, = 0, R-J 2 = R^2 T-j = I which leads 

to: 


TBS=TgO-R 


a2^ 


= T Co 
g a2 

where in this case ' e ^2 same as the emissivlty of saline soil (Table 1). 

-Since a varies directly with frequency, the effect of attenuation by the top layer 
• would be expected to be negligible at 1 .4 GHz in comparison to 13.9 GHz and, 
■especially, 19.35 GHz. Moreover, since h, the depth to brine, varies between 
-about 10 cm near the center of the desert floor and 2.7 m at the margins, the effect 
-of this variation should be reflected in the difference between the measured bright- 
-ness temperature Tg and the brightness temperature of saline soil, T^ Due -to 
the lack of exact information on the moisture, salinity and temperature profile^^'bf 
■'the Great Salt Lake Desert floor, only qualitative comparisons can be made between ■ 
■the measured brightness temperatures and the inferences drawn from the two-layer 
mode! discussed above. 

Table 2 is a summary of vehicle, sensor, and test site parameters pertinent to 
the present study. More detailed information is given in each of the following sub- 
sections. 

3.1 S-193 Observations 

Three different microwave. sensors operating at the same frequency of 13.9 GHz 
ere incorporated in S-193: a radiometer,^ a scatterometer and an altimeter, all sharing 

the same antenna and some of the receiver front-end. In general, for a given pass (or 
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Table 2, Summary of Vehicle, Sonjor, Pass and Site InformaHon. 


S-I94 


S-T93 RAD 


S-173 SCAT 


Vehicle 

’ Skylab 



Skylab 


Skylab 


AIMfudo 

438 km 


• 

438 km 


438 km 


Sensor 

Rodlometer 


Radiometer 

■ 

Scollorometor 

Frequency 

1,4 GHz 


13.9 GHz 


13.9 GHz 


Effectivo An^enno Boomwldth 

15.2° 



2.02° 


1.56° 


Beam Seonning 

None 



Machanlcal 


Meehan Icol 


Pass Information 




• 




Pass Number 

5 

16 

5 


16 

5 

T6 

Poss Dote 

6-5-73 

B-8-73 

6-5-73 


8-8-73 

6-5-73 

8-8-73 

Local Timo 

10:58 a, m. 

9:01 a.m. 

10:58 a.m 


9:01 a.m. 

10:58 a.m. 

9:01 o.m. 

Mode 

nodir looking 

nadir l.ooking 

CTC 


CTC 

CTC 

CTC 

PUch Anglo 

0° 

0° 

0° 


17.0° 

0“ 

17,0° 

Scon Anglo 

N/A 

N/A 

+ 10.5° 


+10.5° 

1 10.5° 

+ 10.5° 

Nadir Anglo Range 

N/A 

N/A 

0°- 10.5° 


l7.0°-20.0° 

0°- 10.5° 

17.0°-20.0“ 

Ground Resolution at 0°-Scon Anglo 

1 17 km diameter 

117 km diameter 

15 ,5 km diometer 

16.2 km X 16,9 km 

11.11 km diameter 

11.6 kmx 12.1 km 

Ground Resolution ot Etsd of Scon 

N/A 

N/A 

15.7 km X 16.0 km 

16,5 km X 17.5 km 

1 1 ,3 km X 1 1 .4 km 

1 1 .8 km X 12.5 km 

Tost SUe Ground Temperoture 

296°K 

298° K 

296°K 


29B'’K 

296°K 

298°K 


General Reference 


NASA (1974) 


Sobtl (1973), 


Sobtl (1973) 


tSMR 
Nimbus 5 
1100 km 

Rodiomeler 

19.35 GHi 

• 

Electrical 


12:00-1:00 p,rt«. 
CTC 
0 ” 

■4 50° 

0°-50° 

25 km diameter 
160 km X 45 km 

Wllheit (1972) 


*$ee section 3.3 

**CTC = Croiss-Track-ConttguQui 



parts thereof) either the altimeter or the RADSCAT (Radiometer“5catterorneter)wos 
operated. For the RADSCAT portion several^modes of operation were available; each 
mode specifies the data taking sequences in terms of pitch and roll angles, polarization 
and sensor. Because of either coverage', incident angle and/or scanning mode considerations, 
only two S“193 passes were judged suitable for the purposes of the present study. These 
are Pass 5 on June 5, 1973 and Pass 16 on August 8, 1973. Detailed discussion of the 
data recorded on Pass 5 will be presented next followed by a brief summary of the 
observations noted from the data of Pass 16. 

Pass 5 was a North~West to South-East descending pass during which S-193 RAD- 
SCAT was operated in a cross-track contiguous mode with the radiometer and 
scatterometer measurements interlaced in time. Other mode parameters include: 
linear polarization in a direction parallel to the spacecraft velocity vector, approximately 
pitch, and roll scan between +10.5° and --10.5° relative to nadir in the cross- 
track direction. In this mode, the signal polarization relative to the ground is 
horizontal. During each scan 12 data points are recorded by each sensor. Figure 5 

shows the radiometer footprints (calculated on the basis of the antenna beamwidth) 

* . ^ 

;enc!osed in the test site frame chosen for this investigation. For reference considerations, 

the scans have been designated by alphabetical letters and the footprints within each 

scan are numbered. The position accuracy of the center of a footprint is estimated to be 

-approximately equal to one half of a footprint diameter. 

In Figure 5, scans a-d, which fall completely outside the desert area, show bright- 
ness temperatures between 270°K and 279°K. As the spacecraft moved towards the south- 
east, lower Tg values were observed over the desert, the lowest being 200°Kat footprints 6h 
and 8h, situated approximately- in the central part of the desert. Furthermore, the 
radiometer shows a drastic change as it crosses boundaries; viz., 3g(264°K) outside 
the desert, 4g(220°K) at the boundary, and 5g (206°K) inside the desert. In scan 
m, as the antenna scanned between 10m and 12m, two different boundaries are crossed; 
lOm is at 260°K, 11m is 205°K (saltflats) and 12rn Is 154°K (lake). Keeping in mind 
the position accuracy of these footprints and the topographic accuracy of the terrain 
map, it is clear that the radiometer response is fairly consistent throughout. Consider, 
for example, footprint lOS. Its relatively low brightness temperature of 198°K leads 
us to suspect that a larger portion of Utah Lake is enclosed in lOS than is shown. 
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S~193 radiomei’er fooi•prin^s, Pass 5, 6/5/73 








For a closer undersl'andlhg of the spatial variation of the brightness 
temperature within the Great Salt Lake Desert area. Figure 6 shows a contour map 
generated on the basis of the data shown in Figure 5.. Based on our earlier discussion 
of the desert hydrology, we propose that the contour map represents a gross "mapping" 
of the depth to the brine layer beneath the surface. 

Whereas the microwave emissivlty is almost insensitive to incident angle 
variation between nadir and 10.5°, the radar scattering coefficient experiences -its 
largest drop in magnitude-in precisely this range of angles. Hence for the purposes of 
this study, only the data within the + 2.5° range around nadir was included 
in the contour map of Figure 7. The radar data confirms the observations made by 
the radiometer in that its measured scattering coefficient changes from values smaller 
than 0 dB for footprints outside the desert to values as high as 17.7 dB over the central 
part of the desert. 

During the early stages of this study, it was suggested that perhaps the response 
observed by S“193 over the Great Sait Lake Desert might be due to spots ;bf5tanding 
water formed by a previous rainfall . No such evidence was. observed-on the Skylab 
photography (S 190A) of this area. The next step was to determine whether or not 
this phenomenori is also evident from data acquired by a) other Skylab passes or 
-b) other satellites carrying microwave sensors. The answer to the latter was provided 
•by ESMR (section 3 .3) while the answer to the earlier was provided by Skylab Pass 16 
•-on August 8, 1973, approximately two months after Pass 5. The subsatellite track. of 
Pass 16 was approximately parallel to that of Pass 5 but was shifted towards the North- 
East by about 52 km. 

During Pass 16 over Utah, S-193 RADSCAT was operated in a cross-track 
contiguous mode similar to Pass 5, except that the pitch angle was about 17° (0° for 
Pass 5). Corresponding to a roll scan of + 10°, the Incident angle varied between 
.about 17° (at 0° roll) and 20° (at 10° roil). Thus, for all practical purposes, the 
incident angle was almost a constant. Due to the larger pitch angle employed 
in Pass 16, the scatterometer and radiometer footprints were slightly larger (than those 
of Pass 5). The Pass 16 data were used to generate footprint and contour maps similar 
to those shown in Figures 5 and 6. Over the test site area, the radiometer recorded 
a high of 273°K and a low of 199°K [Figures 8 and 9 ]. The general shape of 
the contours are similar to those produced from the June 5, 1973 data. Footprints 
over the Salt Lake itself recorded brightness temperatures as low as 146°K. The point 
to'be made here is that indeed the phenoijienon observed in June was also evident in 
August. 
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S-193 BRIGHTNESS TEMPERATURE CONTOURS 
GREAT SALT LAKE DESERT 6/5/73 
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Figure 6, S-193 brighlriess temperature contours based on the data of 











Figure 8. S-193 radiomril'or (oofprini’s. Pass 16, 8/8/73 
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Figure 9. S-193 brighrness temperalure contours based on the data 
of Figure 8, 




3.2 5-194 Obsen/aHons 


Simultaneous with the S-193 measurements discussed in the previous section, 
passive L-band (1 .4 GHz) data were acquired during the same passes over Utah,' The 
radiometer was nadir-looking having a circular footprint approximately 117 km in 
diameter. Ground distance covered between consecutive measurement points shown in 
Figure 10 is approximately 7km, producing about 94 per cent overlap. The measured data 
are plotted in Figure 1 1 os o function of distance from the center of a reference foot- 
print, This plot represents a convolution of the ground brightness temperature (Tg^) 
spatial variation with the antenna power pattern as the latter is swept across the test 
site along the subsatellite track. Before crossing the desert boundary on the western 
side, the S-194 radiometer recorded a constant brightness temperature of 260°K for 
footprints 1 through 9. Between footprints 9 and 26, whose centers are separated by 
approximately 115 km, the brightness temoerature dropped to a iovy of 202°K, after . 
which it remained unchanged over the next 20 km and then it Increased to a value of 
,252^ over the eastern section of the test site. With a 94 per cent overlap between succes- 
-sive footprints, the fast rate of change of the brightness temperature as the antenna beam 
swept across each of the two boundaries of the desert, signifies a much faster rate of 
change of the ground brightness temperature (Tg^) spatial distribution in the direction 
of the satellite track. The combination of S-193 and S-194 observations suggests 
-that the lowest Tg^ value of the desert center at 1 .4 GHz must be much lower than 
202°K, 

Following the same basic procedure as described above, data recorded on 
Pass 16 (August 8, 1973) are shown in Figures 12 and 13, Figure 13 shows a plot 
similar in shape to the plot shown in Figure 11 although having a slightly wider 
range of Tg; west of the desert Tg®=' 270°K, east of the desert Tg®“ 260°K and over 
the desert the lowest Tg recorded was 195°K. The difference between the brightness 
temperatures corresponding to footprints on the Western and eastern sections of the 
test site (10°K) is comparable with the difference observed for Pass 5 (8°K), yet the 
absolute value is higher by about 10°K. Since the teported air temperature for Pass 
16 is only 2°K higher than that of Pass 5, the 10°K increase in absolute value must 
be, for the most part, attributable to differences in the topography due to the ground 
separation of the two passes. The low of 195°K is attributed to partial contributions 
by the Great Salt Lake. 
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S~194 radiometer footprint centers, Pass 5, 6/5/73. Only footprints 
6, 7, 22, 38 and 51 ore shown. 
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S-194 BRIGHTNESS TEMPERATURES 
OVER THE GREAT SALT LAKE DESERT 6/5/73 
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Figure 11 . S-194 brightness temperature as a function of distance from the center 

of footprint 1 on Figure 8. 
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Figure 12. S“194 rad^ome^er foofprin)- cenl-ers, Pass 16 , 8 / 8 / 73 , ' Only fool-prinfs 
1, 18, 37 and 54 are shown. 




S-194' BRIGHTNESS ’ TEAAPERATURES 
OVER THE GREAT SALT LAKE DESERT 8/8/73 



In support of Pass 16, aircraft data was acquired by NASA using the L~band 
portion of the Muiti Frequency Microwave Radiometer (MFMR) on August 10, 1973, 
two days after the test site was traversed by Skyiab . The course flown was a straight 
line from 41,275°N, 113.899'\V to40.843°N, 112.992°W, which crosses the 
desert along a line close to the subsatellite track of Pass 16 [NASA, 1974]. 

Figure 14 is a Skyiab photographic image of the northern part’of the Great 
Salt Lake Desert on which the course flown by MFMR is indicated. Also shown is a 
plot of the brightness temperature recorded by the MFMR instrument which indicates 
Tg values ranging between 285°K and 303°K over the mountains on Ihe two ends of 
the pass and over the New Foundland mountains approximately haif~way in between. 

In comparison, the brightness temperature of the desert portions of the pass are generally 
lower by 150°K or more. The lowest recorded Tg value was 93°K. 

3.3 'Nimbus 5 ESMR Observations 

The Nimbus 5 satellite launched on December 11, 1972 carries an Electrically 
.. Scanning Microwave Radiometer (ESMR). ' . 

The ESMR consists of the following major components: [Wilheit, 1972] 

a) A phased array microwave antenna consisting of 103 
•waveguide elements each having its associated electrical 
phase shifter. The aperture area is 83.3 cm x 85.5 cm. 

The polarization is linear, parallel to the spacecraft 
velocity vector, 

b) A beam steering computer which determines the coil 
current for each of the phase shifters for each beam position. 

c) A microwave receiver with a center frequency of 19.35 GHz 
and an IF bandpass of from 5 to 125 MHz; thus it is sensitive 
to radiation from 19.225 to 19.475 GHz, except for a 10 MHz 
gap in the center of the band. 

The unit is arranged to scan perpendicular to the spacecraft velocity vector from 
50^ to the left to 50° to the right of nadir (in 78 steps) every four seconds. The beam 
width is 1 .4° x 1 .4° near nadir and degrades to 2.2° crosstrack x 1 .4° downtrack at the 
50° extremes. For a nominal orbit of 1100 km altitucfe, the resolution is 25 km x 25 km 
near nadir degrading to 160 km crosstrack x 45 km downtrack at the ends of the scan. 

■ The satellite orbit is polar with noon and midnight equator crossings. 
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MFMR BRIGHTNESS TEMPERATURE VARIATION 
ACROSS THE GREAT SALT LAI^E DESERT. UTAH 8/10/73 



• Figure 14. Skylab photogrcphic Image of the Northern portion of the 
Great Salt Lake Desert on which MFMR coverage on 
August 10, 1973 is indicated, 
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There Is an empirical correction applied to the data for the off nadir beam 
positions, this correction is designed to eliminate the angular variation of bright~ 
ness temperature over the ocean. To minimize the effect of this correction and 
the degradation of resolution, the analysts in this paper was primarily restricted to 
-data from nadir angles of 30° or less. 

A contour map of the ESMR brightness temperatures for the Great Salt Lake area 
is shown in Figure 15 with the out'lTnes of the lake and desert indicated for reference . The 
brightness temperatures over the lake are not as low as those observed over the ocean (less 
than 140°K), because the width of the lake Is comparable to the resolution of the Instru~ 
ment. Thus it is unlikely that the minimum over the desert is as low as it could be. The 
minimum Tg observation over the desert is about 70°K less than the surrounding areas, 
which were up to 280 K. On the following night, the minimum over the Great Sait 
Lake Desert was about the same while the observed brightness temperature over the 
surrounding desert dropped to below 260°K. Similar maps were produced for the area 
on a several per month basis from June, 1973 to December, 1974. The minimum 
brightness temperature observed over the desert is plotted as a function of time for that 
period (Figure 16), For comparison the brightness temperature of a spot (40°45'N x 
114°45‘V/) 60km v/est of Wendover is also plotted to indicate the type of seasonal 
variation that might be expected for the usual terrain outside the desert area. In 
general the brightness temperature for this location varies rather smoothly from a 
maximum of 280°K in July and August down- to a minimum of 240°K for January and 
February and appears to repeat from one year to the next. The minimum brightness 
temperature over the Great Salt Lake Desert followed a similar sort.-of seasonal variation 
with its maximum occurring in July and Augusts However, the minima were in November 
in response to the fall rains during both years. There is a significant difference in the 
level of the brightness temperature for the summers of 1973 and 1974. The minimum 
temperatures observed over the Great Salt Lake Desert in the summer of 1974 were 
20°“30°K higher than those observed during the summer of 1973. The rainfall 
during the summer of 1974 was only 50 per cent of normal while in 1973 It was slightly 
above normal. The average monthly rainfall for eight stations surrounding the basin 
is shown as bar graphs on the bottom of this figure and clearly indicates this difference 
in rainfall for the two summers. The response to the heavier fall rains of November, 

1973 and October, 1974 is Indicated and in particular the lowest brightness temperature 
was observed on 18 November 1973 when more than 1 cm of rain was recorded at the 
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ESMR BRIGHTNESS TEMPERATURE CONTOURS 
GREAT SALT LAKE DESERT 6/5/73 
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Figure 15, ESMR brightness temperature contours, 6/5/73, 
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MONTHLY 

RAINFALL BRIGHTNESS TEMP., KELVINS 

IN cm' 



1973 1974 

Figure 16. Temporal variaHons of ihe minimum recorded ESMR brighf-ness femperafure 
over fhe Greaf Sal}' Lake Deserf (x) compared wiih fhe ESMR brlghfness 
femperature (o) of the reference point outside the desert {indicated in 
Figure 11). 



eight stations. It appears that the radiometer may be responding to a combination 
of surface water resulting from the rain and a rise in the level of the water table. 

There are several other areas of the world which have similar geological 
features; e.g. the salt deserts of Bolivia: the Solar de Coipasa and Solar de 
Uyuni. Figure 17 is a ESMR brightness temperature. contour map of this area, for 
6 June 1973 . Again the salt deserts have a much lower brightness temperature than 
the surrounding areas. These salt deserts are at the southern and lower end (altitude = 
3660 m) of the drainage system of the great central plateau or Altiplano region of 
Bolivia. Thus the topography of these deserts is similar to that represented in Figure 2 
with high mountains (5000-6000 m) on the east and west sides of the plateau, • 

Because of the size of Solar de Uyuni (100 x 120 km) it Is possible to have' ESMR- ' 
resolution elements looking only at the Salar. The lowest brightness temperature 
observed during this pass over the de Uyuni was 165°tC. This is comparable to the 
lowest value observed over the Great Salt Lake Desert after a rain and is comparable 
-to the-value calculated for a saline soil . It is unlikely^ however, that the low value 
.-over the de Uyuni was after a rain since the station at Oruro (200 km north of the 
de Uyuni) reported no rain for June, 1973. Thus in this case, the radiometer is most 
-likely responding to subsurface water. 

‘4.0 CONCLUDING REMARKS 


Comparison of the microwave data acquired by Skylab, Nimbus and MFMR 
with the available information on the hydrology of the Great Salt Lake Desert suggests 
the following conclusions; 

a) Subsurface water exhibits a strong influence on the measured - 
scattering coefficient and brightness temperature, particularly 
over the central part of the desert where the water table is the 
closest to the surface. Although this influence is clearly discernible at 
13,9 GHz and 19,35 GHz, it is most pronounced at 1 .4 GHz; 
over the desert the minimum recorded brightness temperature by 
the L“band channel of the airborne MFMR radiometer was 93"’K 
in comparison to 285°K”303°K over the mountains. 
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ESMR BRIGHTNESS TEMPERATURE CONTOURS 
BOLIVIAN SALT DESERTS 6/6/73 
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Figure 17. ESMR brighfriess J-emperarure con^ours of ^he Bolivian salf deserts, 
6/6/73. 
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b) According to the calculated emissivity of saline soli (with 150 °/oo) 

at 1.4 GHz (Table 1), it is not surprising that the 93°K brightness 
temperature recorded by MFMR is lower than the brightness temperature 
of water at the same frequency and ground temperature. For == 
296°K, = d0.65*^K. The presence of the surface layer 

above the brine increases ^BS to a higher temperature . 

c) The configurations of the brightness temperature contours of the S-193 
radiometer and the Nimbus 5 Electrically Scanning Microwave 
Radiometer are in agreement with the reported variation of the depth 
of the water table below the surface between the central part of the 
desert and the desert margins. 
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APPENDIX E 

PRELIMINARY RESULTS FROM SL-4 (WINTER) MEASUREMENTS 
WITH THE SKYLAB S-193 RADIOMETER/SCATTEROMETER 


The flight occupancy of Skylab was during three periods between May, 1973, 
and February, 1974. The winter occupancy, mission SL’"4, provided the first oppor*- 
tunity to measure significant amounts of snow-covered ground with the 13.9 GHz 
S“193 scatterometer. Unfortunately, a malfunction in the scan mechanism of the S-193 
during September, 1973, resulted in degraded winter observations. During the partial 
repair of the scanning mechanis0i_f .the antenna feed structure was damaged, with a resultant 
large reduction in antenna gain and great increase in side-lobe level. The width of 
the main beam was increased somewhat and its shape altered. This problem made the 
radiometer data useless, since more energy was in the large main sidelobe than in the 
central beam; but the scatterometer could be used because of the smaller sidelobe response 
resulting from the product of receiving and transmitting gains. Because of these problems, 
however, the arrival of the SL4 data at The University of Kansas was delayed until 
March, 1975, and It arrived in a form requiring subsequent processing before it could 
be used by the investigators, v/ho received it in June, 1975. Consequently, this brief 
note only gives an indication of the SL-4 data, and detailed analysis is reserved for a 
later report. 

The damage to the scanning mechanism was repaired by the astronauts in an 
EVA, but it was only possible to repair the cross-track scanner. Consequently, all 
SL-4 measurements were made with Pitch 0^. On January 9, 1974, a further failure 
in the scanning mechanism occured, so that It could scan only to the right, but not to 
the left of the Skylab. Unfortunately, some of the best snow measurements had to be 
made with the system operating in this crippled condition. Most were made in a nominal 
Pitch 0, Roll 0 mode. 

Scatterometer Calibration Corrections 

Much effort was expended by a multi -institution team to modify the processing 
for the SL-4 data from the techniques used for SL-2 and SL-3. Necessarily this in- 
cluded some rather gross estimates, for the antenna pattern had to be reconstructed 
partly from a few measurements made on the ground by the University of Kansas after 
the failure and partly from Inducing a similar failure in the backup antenna — there was 
no assurance the induced failure was exactly the same as the failure in space. The 
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resulting corrections brought the scatterometer data into a relatively correct set of 
values, but further refinement in these values was accomplished on this project by a 
comparison of SL2/3 measurements with SL4 measurements from the same type of 
target area . 

Since the scatter from the oceans is believed to be largely independent of 
season, except perhaps in areas where the temperature of the oceans approaches 
0° C, the ocean seems to be the best candidate for comparing SL2/3 with SL4 
observations. Estimates of the surface winds were provided for our Skylab oceanographic 
studies by personnel at City University of New York under Prof. W. J, Pierson, Jr. 

For a comparison between the missions to be valid, it must be made for the same wind 
speeds. 

Because of uncertainties about the exact wind speeds at most individual points 
on the oceans, caused by the scarcity of ships equipped with adequate measuring 
instruments and making frequent enough observations, single measurements of wind 
speed and scattering coefficient were not compared. Rather, a scattergram was prepared 
relating measured scattering coefficient to wind speed estimates from ship reports. 

Since many observations were made at wind speeds near l8 knots, the values of scattering 
coefficient that best fit the scattering vs. wind“speed relations at 18 knots were determined 
for SL2/3 data and separately for SL~4 data. These determinations were made at 17° 
and 33° angles of incidence, and the differences between the two observations used to 
correct the SL-4 scattering measurements. These observations are reported in Table 
E-1, and the difference between SL”4 observations and SL2/3 observations is plotted 
in Figure E"1 for the VV and HH polarizations as a function of angle. 

TABLE E - 1 

COMPARISON OF SL2/3 AND SL4 RESPONSE OVER OCEAN 
WITH 18 KNOT WINDS 


Polarization 

Angle 

5L2/3 

SL4 

SL4 - SL2/3 


(°) 

a° in dB 

a° in dB 

cr° in dB 

VV 

17 

0.4 

6.7 

.6.3 

HH 

17 

0.4 

5.1 

4.7 

VV 

33 

-9.6 

-2.6 

7.0 

HH 

33 

-11.0 

-5,6 

5.4 





Near vertical the situation was somewhat more complicated. The uncer- 
tain change in the pattern of the antenna, combined with the extreme sensitivity of 
the scattering coefficient to angle within a degree or so of vertical over ocean, make 
the use of vertical incidence data for determining corrections questionable. Unfortunately, 
the CTC mode, which would have provided ample points at angles out to 10° or so, 
was not used over the ocean during SL-4. Consequently, an effort was made to find a 
suitable land surface that had been observed on both winter and summer missions, that 
had very little vegetative cover during the summer and no snov/ in winter so Its winter 
response should be similar to the summer response, and that was observed within a 
few degrees of vertical . This turned out to be more difficult than at first thought be- 
cause the CTC pitch 0° summer passes were mostly In areas covered by snow In winter, 
or in areas of heavy summer vegetation. 

Fortunately, suitable angles were found in- two ITC passes over dry rangeland 
In the Texas panhandle during SL2, and a CTC pass was made over The same area 
during SL4. On day 156 (5 June 1973) an ITC HH-polarization pass went from Utah 
across parts of Colorado, New Mexico and Texas, On day 162 (1 1 June 1973) 
crossed Colorado, Oklahoma, and Texas, On day 11 (11 January 1974) during SL4 
a CTC Pitch 0° W"poIarization pass was made over Texas, Oklahoma, Kansas and 
Missouri, All of these passes covered about the same area of rangeland in the Texas 
panhandle, and the values of scattering coefficient for the almost-common area were 
separated out. Mean values and standard deviations were calculated for these data, and 
the SDs were small enough to warrant using the data for this comparison. Although both 
VV and HH passes were used for SL2, only a VV pass was available for SL4. During 
SL2/3, however, the calibration for VV and HH was the same, so that both of these 
could be used to cross-check the VV from SL-4; this is possible because for angles 
within 10° or so the only difference between VV and HH is the azimuth direction of 
the electric vector and this is not important over such a large area. The resulting 
difference was found to be about 5.6 dB at a nominal 2,3°. This is plotted on 
Figure 1 and shown to be consistent with the trend from the 17° and 33° oceanic 
data, so the corrections Indicated in the figure have been applied to the SL-4 data. 

Actual values used for correction are shown in Table E“2. 
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Correction in dB - Subtract this amount 
from SL4 Data Provided by NASA 


SL4-SL2/3 VV 



Figure E“1 . Corrections To Be Made To SL^ Scotterometer Data, Points 
At ] 7 ^ and 33° Based on 17~*Knot Winds. Point At 2.3° 
Based on Texas Parihandie Fiangeland. 




TABLE E-2 



Corrections Applied to SL-4 Scattering 

Coefficients 

Angle 

W Correction 

HH Correction 

o 

o 

o 

- 5.8 dB 

-4.2 dB 

10 -20° 

- 6.2 dB 

-4.6 dB 


The few values for larger angles are obfained direcfly from the curve. 

Some Preiiminory SL~4 Results and Comparisons with SL2/3 

SL“4 observations have been reduced and analysis and interpretation are 
underway. In this section, we present three examples of data plots that can be 
produced from the overland SL""4 data, along with some comparisons with SL2/3 
data. No attempt should be made to generalize from these curves, for generalization 
will require much more study of the complete body of SL“4 scattering measurements 
and the related ground truth. In particular, conclusions should not be drawn on 
the basis of this limited sample os to the effect of snow. An attempt is being made 
to determine this effect In analysis of the complete data set, but it is not easy to do 
because of the difficulty in finding measurements made under conditions comparable 
except for the presence or absence of snow. 

Figure E“2 shows examples of SL~4 returns from the same Iand“use category, 
Agriculture-Forest (mixture within a resolution cell)* Quite clearly, the nature of 
the forests in the three cases will be different, although agricultural terrain in the winter 
should be relatively free of vegetation in all areas, with the possible exception of 
Louisiana-Mississippi . Perhaps’ this explains why the Pass 87 curve is above the others 
— it represents scatter from standing crops and forests that have significant amounts of 
moisture in the branches and trunks — and in some places the leaves may not have 
fallen. Most of the forest on the Texas-West Virginia pass is far enough north so that 
the deciduous trees will have shed their leaves and the trunks and branches will be 
relatively dry, and the same is probably true for the pass from the Texas panhandle to 
Iowa; actually the snov/-covered part of this pass reported here is from southern Kansas 
to Iowa across part of northern Missouri. The lov/er value of scattering for this pass is 
too close to that for the non-snow pass on day 14 for the difference to be significant 
in terms of snow — it could very well have been just this different without snow at all. 
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Agricultu.re/Forest ( AF) 
Polarization: VV 

Snow, Pass 81, Day 11, 
Texas / Iowa 

v™ Non-snow, Pass 83, Day 14, 
Texas/W. Virginia 
Non-snow, Pass 87, Day 21 , 
Louisiana/Mississippi 




4 6 8 1C 

Angle of Incidence (Degrees) 


Figure E~2. Example of a Single-Category Set of Summer and Winter Measurements Including One SnowrCovered Area 
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Non-snow Winter vs. Summer 
Polarization: VV 

Non-snow, Pass 87, Day 21, 

SL4 Arkansas /Gulf of Mexico 

Non-snow, Pass 83, Day 14, 
Texas/Pennsylvania 

SL3 m— Day 253, CTC, Colorado/Nebraska 
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Figure E"3. Examples of Measurements Without Snow in Winter Compared with Summer Measurements in Different Areas, 



°{dB) 



Angle of incidence (Degrees) 


Figure E“4. Example Comparing Both Snow-Covered and Snow“Free Winter Measurements with Summer Measurements. The 
Wyoming/Nebraska and Arizona/Nebraska Measurements Include Similar Areas, But the Latter Also Includes 
Mountains. 



Figure E“4 compares one summer pass with 2 snow and 2 non-snow pass 
sections from SL-4. Note that in this case the summer returns are all above the 
winter returns, which is the opposite of what happened in Figure E-3. The 
difference between snow and non-snow is too small to be significant. 

The primary conclusion one can draw from such a presentation is that the 
SL4 data appear to be in general consistant with SL2/3 data after we make the cor- 
rections indicated earlier. Firmer conclusions must await more detailed analysis 
which is under way. 
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ABSTRACT 

Certain correlations have been noted between microwave backscatter response 
at approximately 33° from SKYLAB's S“193 Scatterometer (operating at 13.9 GHz) 
and the configuration of ground targets in Brazil as discerned from coarse"^cale maps. 
In spite of the limited nature of available ground truth, these correlations have been 
sufficient to permit the production of image~Iike displays which bear a marked re~ 
semblance to known terrain features in several instances. These terrain features have 
been identified by maps because S-190A and S“190B photographic imagery was almost 
useless (the areas being obscured by clouds) . Among the features identified are large, 
nearly homogeneous portions of rainforest including one section of the Amazon basin 
and its attendant riverine formations. 
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Inlroduction 


Backscaftering coefficient measurements over Brazil using the scatterometer on 
the Skylab manned spacecraft were used to create image~iike displays. These image"" 
like displays show a distinct correlation between the backscatter response and the major 
terrain features as discerned from maps. 

The scatterometer on Skylab, designated as S-193 was an interrupted“CW device 
with a mechanically scanning antenna that had an effective beamwidth of 1 ,45°, The - 
scatterometer, operating at 13.9 GHz, shared the antenna and some of the rf hardware 
with a radiometer and an altimeter. The scan of the antenna was predesigned to one 
of four possible sequences which were chosen by astronaut control. Data from only 
one mode (cross track contiguous with a pitch offset of 29°) were used in preparing the 
images as reported here. At this angle the individual cells are ellipses with 13 x 15 km 
minor and major axes. 

The large areas of homogeneous terrain found in South America (particularly 
the Brazilian region) made them attractive targets for producing the Image~like 
display. The ability of the scatterometer to map, with a gross resolution, the surface 
could only be assessed by maps, because the spacecraft imagery obtained was rendered 
almost useless due to intervening clouds. The assessment of the ability of the Skylab 
scatterometer to respond to major terrain feature variations should be tempered by the 
fact that the regions in question have not been authoritatively and exhaustively mapped. 

The Images were created by interpolating the backscattering coefficient data 

ascribed to the middle point of any target to a finer grid. The interpolated values on the 

finer grid were quantized at 0.25dB to one bit and displayed on both a black-and~white 

video monitor and a color video monitor (through a densityto-hue converter). The IDECS*, 

1 

a near reai~time, analog/digital image processing device was used to produce the images. 
Selection of ground truth 

More than 250 frames of imagery which were located over Brazil were received from 
NASA at the KU Remote Sensing Lab. All of these were derived from the Skylab 3 and 
Skylab 4 missions, with the immediate problem that no real time relationship could be 
established between imagery arid scatterometer data. Further, examination of the imagery 
revealed that its utility was impaired by excessive cloud cover (some tracks were completely 
covered for hundreds of kilometers), as well as the exposure and filtering problems experi- 
enced periodically during the missions. Due to this combination of factors, it was decided 

* IDECS stand for jmage Discrimination Enhancement and Combination System. 
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that this imagery should be relegated to a subordinate role in making determinations as 
to the nature of ground targets. 

The potential ground truth base consisted of maps. As much of the territory 
examined was located in the interior, relatively little of it had been mapped at large 
scale. Furthermore, the limited time frame in which the study was to be completed 
made international inquiries, with possible weeks of turn-around time somewhat im- 
practical. Hence, the decision was made to attempt correlations with the coarse 
scale vegetation and physiographic maps already available (often at 1:25,000,000 or 
smaller) supported occasionally by Operational Navigation Charts at 1:1,000,000 scale. 
Selection of the data base 

In the course of the Skylab 2 mission, several passes were made over Brazil 
using the S-193 Radiometer/Scatterometer in Cross-Track-Contiguous, Pitch 29^, mode. 

In this mode, the instrument takes twelve consecutive measurements as It scans from 
side-to-side, each resolution cell or footprint being approximately 153 square kilomet- 
ers in area, with the indivicual scan tracks being separated by approximately 15 kilomet- 
ers (Figure 1). 

The first pass. Pass 7 (Day 161 - 10 June 1973), began off the northeast coast of 
Brazil, near the mouth of the Amazon (Figure 2), and terminated off the east coast near 
Salvador. The second. Pass 8 (day 162 - 11 June 1973) began in northeastern Columbia 
and terminated in southern Brazil, southwest of Brasilia. The third and final pass. Pass 11 
(Day 165 - 14 June 1973) began on the border separating Brazil and Peru, bisecting the 
extreme western tip of the former and terminating in Bolivia. The extraordinary length 
of these passes enhanced their potential value. 

Results of examination of the S-193 Scatterometer data 

The data was visually examined to determine if trends existed which might be 
relevant to the goals of cur study. It became apparent that large segments of the data 
in Pass 7 reflected such trends. In several instances, the scattering coefficients displayed 
remarkable consistency for numerous scans along the pass, giving way to abrupt changes on 
occasion, which suggested changes in surface configuration. In a similar fashion. Pass 8 
exhibited large, nearly homogeneous segments of data with distinctive dips in the value 
of the scattering coefficient at both polarizations, and at several points along the pass. 
Pass 11, when examined, showed traits similar to the other two. 

Quantification of the analysis was carried out by calculating the means, standard 
deviations and upper and lower bounds of the scattering coefficients of selected samples 
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from the data. These samples were taken from the center of several successive scans. 

The centroids of successive sample points, lying with area., identified as homogeneous 
by the Brazilian Forestry Atlas, were positioned approximately Fifteen kilometers opart. 

As expected, these calculations revealed a high degree of homogeneity which supported 
the conclusions arrived at during visual examination and to\a large extent, supported 
vegetation boundaries on which the sampling was based. Illustrating the trends 
mentioned here are the samples of the traces of backscatter from the center points 
of the scans (Figures 3a and 3b). It should be noted that the trends illustrated in 
these samples contrast to the data over North America, where much more variability 
was found. 

Creation of the Images 

The cross"1-rack-contiguous, pitch 29°, mode provides 12 samples (almost 
equally spaced) in a cross-track dimension of 130 kilometers. The spacing between 
the points in the olong-track direction is roughly 16 kilometers. The bockscattering 
coefficient is ascribed to the middle point of the target; adjocent targets have some over- 
lap. The long passes over Brazil were broken up into 48 scan segments for interpolation 
ond subsequent display. The Interpolation to a finer grid was performed by an inverse 
weighting of the distance from a target point with known back scattering coefficient. 

The interpolated values were then quantized into six-bit words providing 0 q 25 dB per 
bit from “12 dB to -4 dB. The range of the bockscattering coefficient occasionally 
was below the lower threshold, and for such values the disploy was equivalent to the 
background. The ocean surface is usually below -12 dB at 33° and consequently appeared 
as the background tone. The calibration wedge was prepared by assigning values from 
0 to 32 bits in one bit intervals. The files and the calibration wedge were displayed 
on a block-and-white, and on a color display through a density-to-hue converter which 
is part of the ID ECS system. The display was produced under computer control. 
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Correlations between the S~193 Passes ond ‘'Image~like‘' Displays, Pass If Day 161 

This pass began over the ocean (off the Ilha de Mara jo). Since the backscatter 
data were quantified so that scattering coefficients weaker than -12 dB register as black, 
the result is that the first portion of the pass segment is not visible in the image of 
segment 1 (Figure 4), When Figure A is compared to the map of generalized vegetation 
boundaries (Figure 1), it is clear that the general similarity of the predominant vegetation 
was recognized by the sensor. The second segment of the pass continues this general 
pattern in its initial portion, with the image (Figure 5) becoming abruptly irregular in the 
latter portion. At least one vegetation source (Hueck and Siebert, 1972) indicates 

•j 

this area as one containing grasslands with a high water table. This moisture variability 
may have been detected by the sensor, since correlations have been noted elsewhere 
between soil moisture and backscatter coefficients (Uiaby et a!,, 1974). The final full 
segment (of 48 scans) passed through a region generally considered to be dominated by 
Caatinga thornbush vegetation, A high degree of homogeneity was indicated, although 
a lovy'er backscatter Is indicated in the early part of the segment. This lower value probably 
represents an extension of the phenomena seen in the previous segment, or the relative 
openness of the vegetation as opposed to the more closed canopy of the forest areas. Near 
the coast, a belt of coastal forest and semi~dry woodland (transitioned between the former 
and the Caatinga) appeared in the image as a generally lighter area (stronger return). 

Little difference was found between polarizations. 

Pass 8, Day 161 

As in Figure 1, this pass began over Colombia, and passed into Brazil over 
the Rio Icana - Rio (Japes drainage system. The first forty-eight scan segments of the 
pass were over a nearly homogeneous region. The sensor has clearly detected this, as 
the image display for this region exhibited an almost monotonous continuity. 

Not so monotonous is the second pass segment which centers over the Amazon 
Basin proper. Figure 6 shows the image display for this area and may represent the most 
positive correlation obtained within the experiment. The subject area lies almost directly 
to the west of Manaus. In this area, the Amazon undergoes a series of meanders, with 
numerous areas In which the forest opens briefly with lakes of various proportions being 
indicated at intervals. On some vegetation maps, the area directly adjacent to the 
river itself is indicated as a "svv'ampy" rainforest. In any case, some function of the 
riverine system has clearly been "mapped" by the sensor. At two points in the image 
sharp dips in the'backscatter value are noticeable; these areas are dominated by the 
aforementioned lake formations. 
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As in the case of the first segment, the third segment produced an extremely 
homogeneous image display. This was continued into the first half of the fourth and 
last full segment. This area corresponds to a region of rain forest intruded by savanna. 
The last half of the image corresponds to an area which is largely dominated by savanna, 
A lack of homogeneity is indicated In the image. Since navigation charts show abrupt 
changes of elevation in the area, ihe broken nature of terrain and attendant changes in 
vegetation could be expected to result in the variations of backscatter seen on the 
image. 

Conclusion 

It has been demonstrated that, even with the gross resolution of the Skylab 
scatterometer, image-like displays can be produced which clearly show the ability 
of. the scatterometer to recognize variations in major terrain features in areas in 
Brazil where large homogeneous regions exist. The correlation between the back- 
scattering coefficient variations and major terrain type boundaries is high in many 
instances. In other instances, no significant terrain type variations can be found in the 
maps. This implies that the microwave "sees" things that are not necessarily represented 
as variations on a physiographic map, or, that the physiographic maps of the areas 
under study are somewhat outdated and/or erroneous — a not surprising result in view 
of the difficulty of mapping this area by ground or photographic means. Unfortunately, 
the Project RADAM vegetation maps recently prepared in Brazil were unavailable for 
comparison. 
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